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ABSTRACT 


Experimental  and  analytical  studies  are  made  of  the  elastic  and 
elastic-plastic  response  of  a  three  degree  of  freedom  dynamical  system 
composed  of  three  rigid  and  compact  masses  connected  by  slender  columns, 
the  bottom  (heaviest)  mass  being  suspended  as  a  ballistic  pendulum  and 
subjected  to  an  initial  velocity  step.  Such  a  system  is  related  to  the 
behavior  of  machinery  mounted  in  a  ship  subjected  to  an  underwater  ex¬ 
plosion. 

The  analytical  work  was  performed  on  a  high  speed  digital  computer 
and  employed  mathematical  description  of  the  structural  members  in  which 
force-deflection  relations  depended  upon  previous  plastic  strain  history. 

Comparisons  are  made  between  experimental  and  analytical  responses 

* 

on  the  basis  of  the  velocity-t ime  history  of  the  two  smaller  masses. 

Good  agreement  is  shown. 
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Introduce  ion 


The  present  interest  in  the  response  of  structures  subjected  to 
rapidly  applied  dynamic  loads  has  produced  an  Increasing  amount  of 
literature  in  the  area  of  dynamic  loading.  The  literature  covers  the 
behavior  of  structural  metals,  the  response  of  sirple  structures,  and  the 


the  U,  C, 


Naval  Postgraduate  School,  ] xntersy,  California,  dealt  with  the  behavior 


~‘Z  a  sirple  plane  structure  1  -  AiizK  Lae  1  orlzcnl^l  (roof)  'amber  *’£s  of 


great  stiffness  compared  tr  'a  cc1'  ...s  a  sc  i»,  which  the  gicace°t  part 


of  the  rrass  and  weight  was  in  the  tcci.  The  loading  vas  an  impact  in 


the  plane  of  the  structure  and  at  the  rocf  level.  This  systen  could  be 
adequately  approximated  by  a  stngle  iegree  of  freedom  model  mil  a  simple 
theoretical  analysis,  based  essentially  on  energy  considerations,  ras 
able  to  m:,ke  good  predictions  r  '  experimental  behavior. 

In  the  case  of  a  similar  mu / 1 i - s  tory  structure,  the  energy  is  divided 
between  the  various  roof  and  floor  -asses  and  it  is  not  possible  to  use 
energy  considerations  alone  to  obtain  «  theoretical  prediction  of  dynamic 
behavior.  Thus,  analysis  of  a  system  above  a  single  degree  of  freedom 
requires  evaluation  of  the  ji *x  ‘  u  is  equations  of  notion  of  the  system.  If 
the  structure  passes  into  the  plastic  region,  a  suitable  assumption  must 
be  made  to  compute  the  behavior  in  this  region,  fhe  numerical  manipula¬ 
tions  required  in  solving  such  a  problem  become  enormous  and  the  use  of  a 
high  speed  digital  computer  becomes  necessary. 

This  thesis  deals  with  a  small  two- story  structure  mounted  on  a 
large  movable  mass.  The  entire  system  was  reduced  to  a  three  degree,  of 

^Numbers  in |  ]  refer  to  bibliography. 
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freedom  system.  This  particular  system  can  be  likened  to  a  piece  of 
equipment  mounted  to  a  shipvs  hull.  The  ship  is  then  subjected  to  an 
underwater  explosion  and  the  system's  response  is  studied.  Reference  [V] 
gives  examples  of  equipment  and  machinery  damage  and  response,  wherein 
illustrations  of  actual  shipboard  foundation  failures  are  presented  along 
with  a  typical  velocity  history  for  a  mild  shock.  The  ability  to  pre¬ 
dict  the  response  of  such  a  system  would  aid  in  the  original  design  of 
such  structures, 

A  mathematical  model  was  evolved  for  the  above  system  and  a  compari¬ 
son  made  between  the  response  predicted  by  the  mathematical  model  and  a 
test  structure. 

It  is  felt  necessary  to  mention  the  fact  that  the  experimental  por¬ 
tion  of  the  project  was  carried  out  prior  to  the  evaluation  of  the  mathe¬ 
matical  model  and  certain  complications  arose  in  the  mathematical  evolu¬ 
tion  that  would  have  been  eliminated  had  the  sequence  been  reversed.  The 
values  obtained  from  the  mathematical  model  will  be  referred  to  as  the 
"theoretical  results"  throughout;  however,  these  results  cannot  be 
truly  classified  as  "theoretical"  since  a  large  amount  of  experimentation 
with  the  mathematical  model  was  necessary  before  a  suitable  comparison 
could  be  made. 

The  experimental  portion  of  the  thesis  was  conducted  entirely  at 
the  David  Taylor  Model  Basin,  Carderock,  Md.,  during  the  period  30  June  to 
19  July  1961.  The  theoretical  investigation  and  reduction  of  data  was 
carried  out  at  the  U.  S.  Naval  Postgraduate  School,  Monterey,  California, 
during  the  period  1  January  1962  to  1  April  1962,  utilizing  the  Control 
Data  Corporation  1604  computer. 
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The  experimental  investigation  will  be  discussed  before  the  theo¬ 
retical  considerations  due  to  the  fact  that  this  was  the  order  in  which 
they  were  actually  conducted*  A  better  understanding  of  the  complications 
introduced  in  the  theoretical  considerations  will  be  obtained  if  the  factors 
responsible  for  the  complications  are  presented  first. 
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2.  Experimental  Investigation 
Test  Structure 

The  test  structure  consisted  of  two  rectangular  hot  rolled  steel  bars 
measuring  2"  x  2"  x  12",  separated  and  supported  by  two  sets  of  1/8"  x  2" 
hot  rolled  steel  columns.  Figure  (la)  shows  a  sketch  of  th  structure  with¬ 
out  the  attached  measuring  equipment.  The  only  variation  made  to  the  test 
structure  was  the  introduction  of  30Q3-H14  aluminum  columns  for  ten  (10) 
of  the  twenty-six  (26)  test  runs.  These  columns  also  measured  1/8”  x  2". 

The  length  of  the  supporting  columns,  prior  to  mounting  was  12”, 

Holes  measuring  3/8"  in  diameter  were  drilled  on  the  centerline  one  inch 
from  each  end  of  the  supporting  columns  to  accommodate  a  3/8"-NC  bolt  for 
clamping  the  supports  to  the  base  and  bars  mentioned  above. 

The  rectangular  bars  and  base  were  drilled  and  tapped  for  3/8,f-NC 
threads  3/4"  deep  at  both  ends,  to  receive  the  securing  bolts.  Holes 
measuring  3/8"  in  diameter  were  also  drilled  in  the  sides  of  the  bars  for 
attachment  of  the  velocity  meters  used  in  the  test.  Small  metal  bases, 
fitted  to  receive  the  accelerometers  used  in  the  test,  were  welded  to  the 
upper  surfaces  of  each  of  the  rectangular  bars. 

The  entire  clamping  of  the  supporting  columns  was  accomplished  by 
the  single  3/8"-NC  bolt  along  with  a  1/8"  x  2M  x  2"  clamping  plate,  see 
figure  (lb).  The  spacing  between  the  base  to  the  center  bar  and  the 
distance  between  the  center  bar  to  the  upper  bar  was  approximately  8". 

The  single  screw  attachment  of  the  columns  to  the  base  and  bars 
proved  to  be  weak  as  far  as  giving  good  fixity  at  these  junctions.  This 
poor  fixity  caused  the  effective  length  of  the  columns  to  be  somewhat  in 
excess  of  the  actual  free  length  of  8M.  Figure  (2)  is  a  picture  of  the 
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fully  inst rumented  test  structure.  Ihe  velocity  meter  coils  are  on 
the  side  seen  on  the  picture.  The  accelerometers  are  on  the  upper  sur¬ 
faces  facing  toward  the  left. 

The  base  of  the  structure  contained  two  2”  x  2"  x  2"  hot  rolled  steel 
blocks  to  be  used  for  attachment  of  the  supporting  columns.  The  base  was 
designed  to  fit  a  mounting  block  already  attached  to  the  ballistic  pendu¬ 
lum  used  in  the  tests. 

Since  the  width  of  the  supporting  columns  was  large  in  comparison  with 
the  thickness,  the  structured  motion  was  limited  to  a  single  plane. 

Ballistic  Pendulum 

The  ballistic  pendulum  was- used  as  the  means  for  producing  the 
impulsive  loading  on  the  structure.  A  picture  of  the  pendulum  is  shown  in 
figure  (3). 

The  ballistic  pendulum  consisted  of  four  major  components:  1.  The 
anvil,  upon  which  the  test  structure  was  mounted.  2.  The  hammer,  which 
transferred  its  energy  to  the  anvil  producing  the  impulsive  load,  3.  The 
hammer  hold-back  and  release  mechanism.  4.  The  anvil  hold-back  and 
catching  system. 

The  hammer  consisted  of  a  solid  steel  cylinder^  8M  in  diameter  and  48" 
in  length.  Four  smaller  cylinders  were  mounted  on  the  sides  to  receive 
the  supporting  wires  from  the  overhead.  These  cylinders  were  4M  in  d  a  eter 
and  6-l/2u  in  length.  The  total  weight  of  the  hammer  was  785.1  lbs. 

The  anvil  was  likewise  8"  in  diameter;  however,  its  length  was  36M. 

The  four  supporting  cylinders  measured  3,5"  in  diameter  and  7M  in  length. 

The  supporting  cables  were  8  ft  in  lengths  1".  The  total  weight  of  the 
anvil  was  598.6  lbs. 
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ha.*  oi  .c-1  j-back  ci nd  r<L«  se  cehanism  consisted  of  an  electric 


hoist  to  pull  the  hammer  hack  co  an  elevated  position,  the  Initial  posi- 
1 1  on  depending  u'p  ;n  th  e  desire*.  The  refer  sn  mechanism  consisted 

of  an  electro-magnetic  letch  which  could  be  operated  from  a  remote  station 
near  the  recording  instruments,. 

The  anvil  hold-back  and  catching  system  cor.sis.ted  of  a  ratchet  pulley 
and  a  tie-back  lire.  The  line  was  attached  to  the  anvil  at  one  end,  passed 
through  the  pulley  and  carried  25  lb,  weight  at  the  free  end.  Upon  im¬ 
pact  the  anvil  was  forced  back  and  elevate’.  The  25  lb,  v;eight  pulled  the 
line  through  the  ratchet  pulley  and  the  ar.vil  was  prevented  from  return¬ 
ing  upon  reaching  its  maximum  elevation. 

The  ratchet  catching  system  was  used  to  prevent  interference  between 
the  anvil  supporting  wires  and  the  test  structure.  The  proximity  of  the 
test  structure  and  the  anvil  supporting  wires,  when  the  anvil  was  in  its 
equilibrium  position,  was  such  that  the  structure  came  into  contact  with 
the  supporting  wires  upon  impact.  To  eliminate  this,  the  anvil  was  pulled 
back  approximately  8-1/2”  horizontally  prior  to  impact. 

Since  the  experimental  testing  was  conducted  prior  to  a  thorough  theo¬ 
retical  investigation,  the  significance  of  this  initial  anvil  position  was 
not  realized,  and  the  exact  initial  horizontal  location  of  the  anvil  was 
not  recorded*  Upon  examination  of  the  theoretical  aspects  it  was  found  that 
this  initial  position  was  quite  significant  and  complicated  the  theoretical 
solution  of  the  problem.  More  will  be  said  concerning  this  later. 

Me  a suring  Equipment 

The  test  structure  was  instrumented  with  three  bar-magnet  velocity 
meters  manufactured  and  calibrated  by  the  David  Taylor  Model  Basin.  One 
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velocity  meter  was  mounted  on  each  of  the  three  concentrated  masses  of 
the  structure.  The  meter  mounted  on  the  anvil  had  a  sensitivity  of  146 
mv/fps,  the  meter  on  the  center  mass  (lower^ rectangular  bar)  had  a  sensi¬ 
tivity  of  167  mv/fps,  and  the  third  meter,  mounted  on  the  upper  mass  (upper 
rectangular  bar),  had  a  sensitivity  of  151  mv/fps.  From  tests  conducted  by 
the  David  Taylor  Model  Basin,  it  was  determined  that  the  sensitivities 
varied  approximately  it  3%  over  a  4  inch  stroke. 

Two  Statham  accelerometers  were  mounted,  one  each,  on  the  upper  and 
center  masses  of  the  test  structure. 

The  output  from  the  velocity  meters  went  to  the  Velocity  Meter  Control 
Unit  (Attenuator),  manufactured  by  the  David  Taylor  Model  Basin*  No  ampli¬ 
fication  was  possible  with  this-- unit,  the  maximum  output  being  one  times 
the  input.  This  unit  also  permitted  a  convenient  switching  arrangement 
for  obtaining  a  calibration  trace  on  the  oscillograph. 

The  output  from  the  Control  Unit  was  picked  up  by  a  Consolidated 
Electrodynamics  Corporation  Recording  Oscillograph.  The  oscillograph 
utilized  galvanometers  of  the  fluid  damped  type  for  recording  the  velocity 
response  on  the  oscillograph  record.  The  galvanometers  had  a  natural 
frequency  of  1000  cps. 

Outputs  from’ the  accelerometers  were  amplified  by  Consolidated 
Engineering  Corporation  3-KC  Carrier  Amplifiers.  The  power  supply  for 
the  above  amplifiers  consisted  of  a  Consolidated  Engineering  Corporation 
Oscillator-Power  Supply.  The  amplified  accelerometer  signals  were  then 
received  by  the  above  mentioned  oscillograph. 

A  General  Radio  Company,  1000  cycle,  Vacuum-Tube  Fork  was  used  to 
produce  a  time  reference  trace  on  the  oscillograph  record*  The  oscillo¬ 
graph  also  produced  a  time  reference;  however,  it  was  felt  that  the  1000 
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cycle  trace  voulci  ptf.vi.de  a  Oood  cb-  ■*!'  a  t\o  «  scillograpb*  6  tir.<d  sea  la  <, 

Since  difficulty  vas  encountered  with  the  interference  between  the 
test  structure  and  r.he  anvil  supporting  wires,  an  extra  galvanometer  was 
utilized  for  recording  an  indication  of  contact  between  the  wires  and  the 
structure.  This  consisted  of  a  simple  battery  circuit,  using  the  point  of 
probable  contact  as  the  switch.  When  wires  and  structure  made  contact,  a 
pulse  was  recorded  on  the  oscillograph  record,  the  duration  of  this  pulse 
gave  the  period  of  contact.  Figure  (4)  shows  a  sketch  of  the  instruments- 
t  ion. 

Appendix  VI  gives  a  tabulation  of  the  various  instrument  serial  and 
model  numbers. 
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Test  Procedure 


The  following  steps  were  followed  for  a  typical  test  run 

1.  The  structure's  columns  (springs)  were  placed  in  position  and  securely 
fastened, 

2.  The  anvil  was  elevated  to  a  position  where  interference  between  the 
structure  and  supporting  wires  would  not  be  encountered, 

3.  The  velocity  meter  magnets  were  adjusted  to  insure  free  passage  as  the 
anvil  moved  away* 

4.  The  hammer  was  elevated  to  a  predetermined  position*  depending  upon 
the  step  input  desired*  using  the  electric  hoist. 

5.  The  recording  oscillograph  was  operated  at  a  paper  speed  of  4  in. /sec* 
and  velocity  meter  and  accelerometer  calibration  traces  were  recorded. 

6.  All  calibration  settings  for  the  various  instruments  were  recorded. 

7.  The  instruments  were  then  set  for  the  predetermined  settings  for  the 
actual  run,  and  all  values  recorded. 

8.  The  recording  oscillograph  was  operated  at  a  paper  speed  of  40  in. /sec. 
and  the  hammer  released. 

9.  The  oscillograph  was  secured  when  the  velocity  meter  magnets  were  free 
from  the  coils  or*  for  the  lower  input  runs,  when  it  was  felt  sufficient 
data  were  recorded. 
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List  of  Symbols 


a 

b 

E 

F. 

1 

8 

h 

H 

I 

k 

k* 

Jt 

L 

M. 

1 

M 

u 

Q 

R 

T 

V 

x, 

i 

x 

u 

y, 
r 


-  initial  horizontal  anvil  displacerrient 

-  width  of  rectangular  supporting  column 

-  modulus  of  elasticity 

-  force  produced  by  flexure  of  columns  (1  -  1,2) 

2 

-  acceleration  of  gravity  (32.2  ft/sec  ) 

-  thickness  of  rectangular  supporting  columns 

-  (subscript)  denotes  horizontal  component 

-  moment  of  inertia  of  the  cross  section 

-  spring  constant  in  elastic  range 

-  spring  constant  in  plastic  range 

-  column  length 

-  length  of  anvil  supporting  cable 

-  mass  (i  s  1*2,3) 

-  fully  plastic  bending  moment 

-  axial  load  on  columns  produced  by  structures  weight 

-  tension  in  hold-back  line 

-  tension  in  supporting  cables 

-  (subscript)  denotes  vertical  component 

-  horizontal  displacement  of 

-  relative  spring  displacement  at  which  the  fully  plastic  bending 
moment  is  realized 

-  vertical  displacement  of  anvil 

-  dynamic  yield  strength 
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3.  Theoretical  Considerations 


Mathematical  Model 

The  mathematical  model  of  the  two-story  structure  can  be  visualized 
as  a  mass  -  spring  system  as  shown  in  figure  (5a) *  This  idealization 
reduces  the  structure  to  a  three  degree  of  freedom  system*  Figure  (5b) 
shows  the  free  body  diagram  for  the  system*  The  symbols  denote  the 
masses  of  the  anvil,  center  mass,  and  upper  mass  respectively,  the 
symbols  denote  the  forces  produced  by  the  flexure  of  the  columns,  and 
the  denote  displacements  of  the  respective  masses* 

Summing  forces  for  each  mass  results  in  the  following  equations  of 
motion, 


•  • 

X,  =  - 

Ft  +  Th 

M, 

0) 

•  # 

X2  = 

B  -.£* 

Ma 

(2) 

•  • 

*5  ~ 

F* 

m3 

(3) 

where  the  force  T  is  the  horizontal  component  of  the  force  in  the 
H 

support ing  cables* 

Letting  "L"  denote  the  length  of  the  anvil  supporting  wire,  "a" 
the  horizontal  distance  that  the  anvil  is  initially  displaced,  '’x1,1,  the 
displacement  of  the  anvil  from  the  initial  position  "a",  and  "y, ",  the 
vertical  displacement  of  the  anvil,  we  have  from  figure  (6a), 

Lf  =  (°-  +  x .  )2  +  -  J/ 1 ) 

where,  ^  —  VL2— 

In  figure  (6b)  "W"  is  the  total  weight  of  the  anvil  plus  structure 
(W  =  656  lbs,).  The  symbol  "T"  denotes  the  tension  in  the  supporting  cables. 
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The  symbol  ”R"  denotes  the  hold-back  wire  reaction  which  immediately 
reduces  to  zero  upon  impact.  From  figure  (6b)  we  may  write, 


Tv  - 

&  5G  -j- 

Tl  ■ 

a  +  x 

Tv  ” 

Co  & 

H 

1 

from  figure  (6a), 

ili  =  X  —  \/l2  -  (a--tXl)Z 

V  —  (<x+  x,)  x, 

£.  -  y, 

» •  /  v  •  •  •  2.  *2. 

v  —  (OL-t-  x,^  X,  +  X,  +■  V. 


(4) 

(5) 


(6) 

(7) 

(8) 


By  substitution  of  equation  (8)  into  equation  (4)  and  then  combining 

with  equation  (5),  we  arrive  at  a  value  for  T  in  terms  of  known  or 

H 

computed  variables.  Substituting  this  value  of  T  back  into  equation 

H 

•  # 

(1)  and  solving  for  X,  produces, 


x,  = 


-F,  656  G  + 

M.  +  (-^t)2  . 


xf  +  3, 2  I 
9 


(9) 


The  three  equations  describing  the  response  of  the  model  are  equa¬ 
tions  (2),  (3),  and  (9). 


Basic  Considerations 

The  test  structure  followed  closely  to  the  portal  frame  used  by 
Kurzenhauser  in  his  tests,  the  major  differences  being  in  the  number  of 
degrees  of  freedom  and  the  method  of  loading.  The  portal  frame  had  a 
fixed  base  and  the  loading  took  place  through  the  upper  mass  of  the  struc¬ 
ture.  The  two-story  structure  used  in  these  tests  was  loaded  through  its 
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base*  which  dtt'cbtd  to  a  *  bT  wss* 

When  deforrrrt beyond  the  elastic  limit  vert  encountered,  it  be¬ 
came  n ecessavy  to  deviate  fron  Hooke’s  Law  and  to  study  the  behavior  in 
the  plastic  region*  Since  the  entire  system  was  reduced  to  a  spring  - 
mass  system  it  became  necessary  to  assume  a  force  versus  deflection  relation¬ 
ship  that  would  satisfy  the  behavior  of  the  structure  as  it  passed  from  the 
elastic  region  Into  the  plastic  region*  Figure  (7)  shows  the  assumed 
force  versus  relative  deflection  response. 

The  structure’s  columns  behaved  elastically  up  to  the  relative  deflect 
tion  corresponding  to  X*  ,  whereupon  the  behavior  departed  from  Hooke’s  Law* 
The  simplest  way  to  represent  the  jcir.pl icated  response  in  the  elastic-plastic 
defomat ion,  is  to  assu  two  linear  spring  constants,  (slope  of  the  lines 
represented  by  the  force  deflection  curve)  one  for  the  elastic  range  and  a 
second  for  the  plastic  range.  The  elastic  spring  constant  is  represented 
by  k  and  the  plastic  spring  constant  by  k*.  Upon  loading,  the  columns 
follow  the  elastic  spring  constant  (k)  until  yielding  of  the  columns  takes 
place,  then  they  follow  the  plastic  spring  constant  (k*).  When  unloaded, 
the  columns  follow  the  Isplaced  elastic  siring  constant  (k),  forming  an 
hysteresis  loop  as  indicated  in  the  figure. 

It  is  assumed  that  the  transition  from  the  elastic  ^ending  moment  to 
the  fully  plastic  bending  moment  in  the  column  cross-sections,  occurs  in¬ 
stantaneously  at  the  deflection  coiresponding  to  Xu.. 

For  rectangular  cross-sections,  the  fully  plastic  bending  moment  is 
defined  by  the  relation,  2. 

M*= 

where  q=-  is  the  yield  point  or  yield  strength  of  the  material  involved, 
b  is  the  width  of  the  colunn,  and  h  is  the  thickness  of  the  column. 
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A  fotce  analysis  of  one  of  the  columns  in  the  elastic  region  will 
permit  evaluation  of  the  elastic  spring  constant.  Referring  to  figure 
(8a),  and  utilizing  the  cantilever  beam  deflection  equations  gives  the 


following  results , 


_x_  (£■)  (4-): 

2  3  E  I 


=  e  r 

24  EI 


Where:  F  =  total  forcing  function 

on  two  supporting 
colunns 

E  s=  modulus  of  elasticity 
I  =  moment  of  inertia 


Considering  the  influence  of  the  compressive  force  Q,  where  Q  denotes  the 
total  compressive  force  on  two  columns  and  is  the  critical  Euler  column 


loading  for  two  columns,  we  have, 

K  = 


F_  _  24  EI 

X  “  p 


Q 

Q 


A  J 


(10) 


where, 


Q.= 


_  2  TT^E  I 


It  has  been  shown  that  the  approximation  used  in  equation  (10,  ,fis 
sufficiently  accurate  almost  up  to  the  critical  value  of  the  load.  M  [VJ 
A  force  analysis  of  one  of  the  supporting  columns,  figure  (8b),  when 


:s  the 

following 

for 

small 

2Ma 

=  0 

2  M* 

-Q.V 

=  2  X 

I 

i  < 

F  - 

-Q-x 
-  J(  x 

4- 

4  Mm. 

X 

(11) 


Equation  (11)  is  the  equation  of  the  straight  line  representing  the 
plastic  region  shown  on  the  force-deflection  diagram  of  figure  (7)*  The 
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spring  constant  k*  for  the  plastic  region  is  therefore  the  slope  of  this 


line, 


K 


a 


(12) 


i 


It  will  be  noted  from  figure  (7)  that*  at  the  deflection  x^,  the 
spring  force  represented  by  equation  (10)  for  the  elastic  region  is  equal 
to  the  force  represented  by  equation  (12)  for  the  plastic  region.  Equat¬ 
ing  these  two  produces  the  following, 


(13) 


The  value  for  can  be  easily  obtained  from  equation  (13). 


Complications  Introduced  by  Experiments  1  Set-up 

The  two  major  problems  that  arose  from  the  experimental  model  were 
the  failure  to  record  the  initial  anvil  displacement  prior  to  impact  and 
the  lack  of  fixity  of  the  supporting  columns. 

Since  the  initial  anvil  displacement  was  not  recorded,  it  became 
necessary  to  determine  this  value  from  the  available  infer  at  lor.  This  was 
accomplished  by  plotting  the  sum  of  the  system’s  total  momentum  versus  time 
for  both  the  experimental  and  theoretical  values.  The  initial  displace¬ 
ment  "a"  used  in  equation  (9)  was  varied  through  a  reasonable  range  until 
agreement  with  the  experimental  points  was  obtained.  Figure  (9)  is  a 
representative  plot.  The  above  mentioned  comparison  between  theoretical 
and  experimental  was  made  for  each  run  where  sufficient  values  of  anvil 
velocity  were  available.  Since  the  length  of  the  velocity  trace  depended 
upon  the  time  the  magnet  departed  the  velocity  meter  coil,  the  values  for 
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the  anvil  velocity  for  the  plastic  runs  were  not  available  beyond  about 
50  to  100  V  (millis  ~onds).  The  value  obtained  for  the  initial  anvil 
displacement  utilizing  the  above  method  was  consistently  0,7  ft.  The 
value  arrived  at  for  this  displacement  should  have  been  fixed  for  all 
runs,  as  Indicated  from  the  plots,  since  the  anvil  was  moved  to  the  same 
position  for  each  run  during  the  actual  tests. 

The  value  (a  =  0,7*)  was  therefore  used  for  all  runs  with  the  ex¬ 
ception  of  two  of  the  low  input  steel  runs,  which  definitely  indicated  a 

* 

lower  value  for  "a11. 

The  lack  of  fixity  of  the  supporting  columns  introduced  a  problem  as 
to  the  value  for  effective  length  in  the  formulae  for  determing  k  and  k*. 
The  effective  length  was  varied  until  phase  agreement  (peaking  at  the  same 
time)  was  achieved.  The  effective  lengths  found  necessary  were  9,0  inches 
for  steel  and  8,75  inches  for  aluminum. 
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4.  Numerical  Solution 


Reduct  ion  of  Equations  for  Computer 

The  first  step  in  reducing  the  equations  to  a  suitable  form  for 
introduction  into  the  computer  was  to  reduce  the  three  second  order 
equations  (2),  (3),  and  (9)  to  six  first  order  equations.  Introducing 
the  following  dummy  variables, 

f  -  X, 

n  -  x* 
y  =  x3 

then, 

i  =  X, 

•  •  * 

n  =  x2 

•  •  • 

*  - 

the  six  first  order  equations  then  become, 


x,  -  * 

xa=  n 

x3=  r 


-  F,  -  < 

r  OL  -t-x, 
y.  ) 

r.  .  f4+ 

U  +  3  (£-  TT) 

M,  +  ( 

(*  +  Zj  ' 
U-  9,  ; 

r-  ■ 

Gz5Jz 

3 

F,  -  F 
Ma 


(14) 

(15) 

(16) 

(17) * 

(18) 


(19) 


*  Note  that  3^—  j/(,  and  Si  are  given  by  algebraic  relations  involving 
Xj  and  ;  see  Equations  (6)  and  (7),  page  12 


17 


Letting  the  symbol  D  denote  the  dependent  variable,  Y  denote  the 


independent  variable,  and  introducing  the  following  equivalent  symbols, 


D,  =  Y 

Di  =  7 

D,  =  X 

D*  =  j 

Os  =  n 

d6  =  i 


X  =  x, 
y2  =  X, 
y3  =  x2 

=  Y 
Ys  =  7 
Yfc  =  Y 


CC  -  o-  -1-  x , 
CK  =  Y, 
YV  =  S, 


produces  the  following  equations  written  in  the  Fortran  language, 

DO)  =  Y(+) 

D(i)  =  YU) 

CXs)  =  YU) 

DU)  =  (-F(l)  -  (Cc/ck)  *  &5fo.  *(/.+(  Y(4)**2  +  YV**-2)/ 
(52. 20+CK)))/OS-  3+656. /  52.20)  *(CC/CK)  **2) 

DU)  =  ( F( I )  —  F(2))/o.7I  I 

0(b)-  FU)  /  0.705" 


Appendix  I  contains  the  Main  Flow  Chart,  the  Force  versus  Relative  De¬ 
flection  Flow  Chart,  and  the  actual  Fortran  program*  The  program  traces 
the  time  history  of  the  displacements,  velocities,  and  accelerations* 

of  the  three  masses  M. ,  M0,  and  M  starting  with  the  initial  conditions 

\  i.  j 

X1  =  X2  88  x3  =  0;  V2  ®  V3  “  °>  Vi  =  STEP- 


*Accelerat ions  are  given  in  the  actual  computer  print-out,  however,  they  were 

not  used  in  the  analysis  and  have  been  trimmed  off  the  sheets  shown  in  Appen¬ 
dix  V.  Also,  the  print  statement  shown  in  Appendix  I  does  not  produce  dis¬ 
placements.  Displacements  were,  however,  printed  in  an  earlier  version  and 
it  would  be  very  easy  to  restore  them. 
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Numerical  Integraf lor  System 

The  numerical  integration  system  used  was  the  Rungt-Kutta  method. 

The  two  major  advantages  of  this  system,  as  stated  by  Milne  jjpHj  ,  are: 

1.  No  special  starting  procedure  is  required. 

2.  The  length  of  the  time  interval  may  be  altered  during  the 
computation. 

The  reasons  for  utilising  the  Runge-Kutta  method  are:  fisst  of  all  the 
fact  that  the  integration  relies  upon  the  last  computed  point,  as  listed 
in  advantage  number  one  above,  secondly,  the  routine  was  already  available 
for  immediate  use  in  the  U.  S„  Naval  Postgraduate  School  computer  center. 

The  system,  as  with  many  numerical  integration  systems,  is  based  on 
the  Taylor  Series.  For  a  single  first  order  differential  equation  y' ~ 
f(  tr  ,y)  the  following  computations  are  completed  to  arrive  at  a  value 
y (  t  +  At), 

K,  =  At  f  ( t ,  )j ) 

K  »  A-tf  (*  +  3  +  i  k,') 

k5  =  A-t-P  ^  +  i  ^z) 

K4  -  At  f  NJ+  k3) 

^+i[K,  +  2Kit2K3  +  Kt] 

The  same  basic  procedure  as  above  can  be  applied  to  a  system  of  first 
order  equations  such  as  the  six  equations  shown  previously. 
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5,  Graphical  C<:r>pa f i of  results 

This  section  shows  a  graphical  comparison  of  the  experimental  results, 
represented  by  the  circled  points,  and  the  theoretical  response  predicted 
by  the  mathematical  model,  represented  by  the  smooth  curves. 

The  first  section  displays  the  results  of  sixteen  runs  using  steel 
columns.  The  first  seven  plots  contain  a  single  theoretical  curve  since 
these  columns  experienced  no  plastic  deformation.  The  next  nine  plots  dis¬ 
play  three  theoretical  curves  corresponding  to  dynamic  yield  strengths  of 
40,000  psi,  50,000  psi,  and  60,000  psi.  The  last  two  plots  repeat  two  of 
the  runs  extending  the  time  cut-off  to  200  MS. 

The  second  section  displays  the  results  of  ten  runs  using  aluminum 
columns.  The  first  three  runs  are  in  the  elastic  region  and  therefore 
display  the  single  theoretical  curve.  The  remaining  seven  runs  display 
two  theoretical  curves  for  dynamic  yield  strengths  of  20,000  psi  and  30,000 
psi.  The  theoretical  curves  shown  on  pages  44  ,  45  ,  and  46  actually 

merge,  within  plotting  error,  from  approximately  90  MS  to  100  MS. 

The  term  "Step"  denotes  the  initial  anvil  velocity  upon  impact  (ft/ 
sec).  It  was  determined  from  the  velocity  meter  mounted  on  the  anvil. 
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Section  One  -  Steel 
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Section  Two  -  Aluminum  Alloy 
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6.  Discussion  of  Results 


1-  A  complete  tabulation  of  the  experimental  and  theoretical  points  shown 
in  the  graphical  comparison  is  not  included  in  this  thesis.  Samples  of 
the  data  are  contained  in  Appendices  (II)  and  (V).  The  Fortran  program  used 
is  found  in  Appendix  (I),  and,  if  desired,  may  be  utilized  for  obtaining 
additional  theoretical  results.  The  experimental  results  show  on  the 
various  plots  are  as  accurate  as  the  values  used  in  plotting. 

2,  The  two  major  experimental  discrepancies,  as  noted  earlier,  were  the 
lack  of  fixity  and  the  failure  to  record  the  initial  anvil  displacement. 
These  two  items  have  been  suitably  discussed  in  previous  sections.;  however, 
an  improvement  to  the  test  structure  is  suggested  for  future  tests.  The 
single  bolting  arrangement  used  appeared  to  be  inadequate,  as  far  as  pro¬ 
viding  a  fixity  approaching  a  cantilever  beam.  The  more  rigid  connection 
used  by  Kurzenhauser  £l  j  in  his  tests  consisted  of  four  5/16"  N.  C.  socket 
head  cap  screws  plus  a  l/4n  clamping  plate.  The  four  screws  were  located 
near  the  edges  of  the  clamping  plate  instead  of  the  centrally  located  bolt 
used  in  this  thesis.  It  is  therefore  felt  that  the  design  of  the  test 
structure  was  poor  in  this  respect,  and  the  use  of  four  holding  bolts 
would  have  improved  the  fixity. 

Nothing  further  will  be  said  concerning  the  failure  to  record  the  in¬ 
itial  anvil  displacement,  except  that  this  was  vital  information  that  should 
have  been  recorded.  In  any  future  structural  tests  utilizing  the  ballistic 
pendulum  herein  described,  the  test  structure  should  receive  the  initial  im¬ 
pact  with  the  anvil  hanging  vertically  (initial  displacement  zero),  thus 
eliminating  this  initial  horizontal  force. 


51 


3.  The  maximum  step  velocity  (initial  anvil  velocity  upon  impact)  ex¬ 
perienced  by  the  test  structure  was  7.83  ft/sec.  The  ballistic  pendulum 
was  capable  of  producing  approximately  14.0  ft/sec;  however,  due  to  the 
interference  between  the  test  structure  and  the  anvil  supporting  cables, 
the  maximum  of  7.83  ft/sec.  was  the  highest  input  before  interference  was 
experienced. 

4.  The  limiting  value  of  time  used  for  the  graphical  comparison  of  results 
was  100  MS.  The  reason  for  cutting  off  the  comparison  at  this  point  was 
due  to  the  fact  that  in  reconstructing  the  initial  anvil  displacement 

(a  =  0.7'),  the  total  momentum  vs.  time  plot  was  used.  This  plot  was 
limited  by  the  velocity  record  of  the  anvil.  The  velocity  meter  magnet 
for  the  anvil  departed  the  coil  after  approximately  100  MS,  up  to  a  step  of 
3.86  ft/sec.  Above  this  value  for  initial  velocity,  the  anvil  velocity 
trace  was  even  shorter,  averaging  about  50  MS.  Since  the  total  momentum 
vs.  time  plot  was  limited  to  approximately  100  MS,  it  was  felt  that  use  of 
the  predicted  anvil  displacement  beyond  100  MS  would  not  be  warrented. 

A  second  reason  for  using  the  100  MS  limit  was  that  the  two  velocity 
meters  mounted  on  the  structural  masses  produced  faulty  traces  beyond  70  MS 
for  the  steel  columns,  at  the  higher  inputs.  These  faulty  traces  were  like¬ 
wise  caused  by  the  velocity  meter  magnets  departing  their  respective  coils. 
The  purpose  of  mounting  the  accelerometers  on  the  test  structure  was  to 
allow  for  a  continuation  of  the  velocity  plots  beyond  70  MS  by  integration  of 
the  acceleration  vs.  time  traces.  It  was  felt  that  the  increased  error  due 
to  reduction  and  integration  of  the  acceleration  traces  would  be  held  to  a 
minimum  if  the  number  of  points  added  by  this  method  were  small.  Three 
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points  were  added,  producing  a  total  velocity  vs.  time  trace  of  100  MS 
duration.  Discrepancies  of  the  order  i  5%  were  noted  between  velocities 
obtained  from  the  velocity  traces  directly  and  by  integration  of  the 
acceleration  traces,  over  intervals  x^here  the  velocity  traces  were  consider¬ 
ed  valid. 

5.  The  graphical  comparison  of  results,  for  steel  columns,  found  on  pages 
22  through  37  ,  consists  of  approximately  seven  elastic  runs  and  nine 
plastic  runs.  The  experimental  points  for  the  elastic  runs  are  compared 
with  a  single  theoretical  plot  since  the  theoretical  lesults  within  this 
range,  are  not  affected  by  changes  in  the  dynamic  yield  strength.  The  plot 
for  step  -  2.87  ft/s ec  is  a  borderline  case.  A  slight  variation  between 
theoretical  values  forcr^  40,000  psi  and  cr  ^  60,000  psi  was  noted,  how¬ 
ever  the  variation  was  quite  small  (0.3  ft/sec  maximum)  and  therefore  only 
the  theoretical  trace  for  cr =  40,000  psi  is  shown.  Starting  with  the  value 
step  =  3.86  ft/sec,  three  theoretical  curves  are  shown  for  the  remainder 
of  the  runs.  The  curves  represent  dynamic  yield  strengths  of  40,000  psi, 
50,000  psi,  and  60,000  psi.  The  general  trend  in  the  curves  indicates  the 
increased  dynamic  yield  strength  with  increased  initial  velocity.  The  runs 
for  step  =  5.02  ft/sec  and  7.23  ft/sec  do  not  fit  in  x^ith  the  remainder  of 
the  plots.  The  dynamic  yield  strength  for  these  two  runs  shows  a  decrease, 
xdiereas,  an  increase  should  be  indicated  to  follow  the  sequence  of  values 
previous  to  and  following  these  respective  runs.  The  cause  of  the  dis¬ 
agreement  could  possibly  be  due  to  an  error  in  the  experimental  values  for 
initial  anvil  velocity.  The  comparison  indicates  a  lower  velocity  input 
than  the  values  used.  A  change  of  0.02"  in  measurement  from  the  oscillo¬ 
graph  record  produces  a  change  of  0.19  ft/sec.  Also  a  variation  in  support¬ 
ing  columns  from  one  test  structure  to  another  could  introduce  an  appreciable 


error.  It  can  be  seen  in  sore  of  the  higher  input  plots,  for  steel,  that 
the  experimental  points  tend  to  fall  v>ff  more  rapidly  than  would  be  indicat¬ 
ed  for  continuity  of  plot.  As  was  noted  earlier,  the  last  three  points 
for  the  steel  plastic  runs  were  cor, puled  using  the  integrated  acceleration 
vs.  time  traces*  The  added  error  caused  by  this  procedure  may  account  for 
this  discrepancy. 

Pages  38  and  39  show  an  extension  of  two  of  the  previous  runs  out  to 
200  M^.  The  purpose  of  these  two  plots  was  to  give  an  indication  of  the 
agreement  beyond  the  100  MS  limit.  As  can  be  seen  from  the  plots,  the 
agreement  is  very  good.  This  further  indicates  the  validity  of  the  computed 
value  for  initial  anvil  displacement, 

6.  The  graphical  comparisons  using  aluminum  columns  are  likewise  limited 
to  100  MS.  The  first  three  plots  show  the  single  theoretical  trace  ifhereas 
the  remaining  seven  plots  show  two  theoretical  traces,  one  trace  representing 
CT  -  20,000  psi  and  a  second  trace  representing  cr  »  30,000  psi.  For 

the  three  runs  represented  by  step  -  2.29,  2  49,  and  2.82  ft/sec,  the 
theoretical  values  for  sr  -  20,000  psi  and  30,000  psi  merge  within 

plotting  errot  between  90  and  100  MS. 

The  experimental  values  for  the  velocity-time  plot  for  the  aluminum 
columns  were  taken  entirely  from  the  velocity  trace  since  the  magnets  re¬ 
mained  in  the  coils  somewhat  longer  than  for  the  steel  columns.  For  the 
two  highest  input  velocities  the  last  three  points  were  checked  by  integrat¬ 
ing  the  acceleration  vs.  time  traces. 

7.  Massard  and  Collins  [2j  ,  in  tests  of  structural  metals  under  rapid 
loading,  demonstrated  the  time  sensitive  behavior  known  as  "delayed  yield" 
for  structural  steel  and  aluminum.  It  vras  found  that  metals  that  display 
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a  definite  upper  and  lower  yield  point,  under  slow  testing  procedures, 
are  quite  time  sensitive  and  display  the  delayed  yield  phenomenon.  For 
metals  such  as  aluminum,  where  straining  under  slow  tests  is  not  dis¬ 
continuous,  the  delayed  yield  effect  is  not  as  prominent.  From  the 
graphical  presentation  of  data,  the  aluminum  plots  would  tend  to  substan¬ 
tiate  this  conclusion,  since  the  apparent  increase  in  dynamic  yield 
strength  is  somewhat  smaller  relative  to  that  demonstrated  for  steel. 

8.  An  improvement  to  the  testing  procedure  used  would  be  to  replace  the 
supporting  columns  before  each  run.  Obviously,  this  is  a  requirement  for 
the  plastic  runs;  however,  in  several  instances,  the  same  set  of  supporting 
columns  was  used  for  one  very  low  input  elastic  run  and  then  a  higher 
velocity  input  plastic  run.  This  procedure  only  involved  about  three  of 
the  runs  and  was  done  to  save  time. 

9.  As  mentioned  earlier,  the  integration  routine  utilized  in  this  thesis 

was  the  Runge-Kutta  method.  The  time  interval  used  was  1/10  MS.  This  was 
arrived  at  by  comparison  with  the  results  for  an  exact  solution  for  a 
similar  system  contained  in  reference  .  The  only  difference  between 

the  two  systems  was  that  the  exact  solution  was  based  upon  a  system  whose 
base  was  rigid.  The  system  utilized  in  this  thesis  was  easily  reduced  to 
the  above  system  merely  by  making  very  large.  Comparison  of  the  two 
systems  in  the  elastic  range  (exact  solution  was  valid  for  elastic  response 
only)  produced  answers  that  were  in  agreement  to  about  three  significant 
figures.  This  was  considered  accurate  enough  for  plotting. 

Since  the  time  interval  used  in  the  computer  program  was  1/10  MS, 
theoretical  results  could  have  been  obtained  for  every  1/10  MS;  however, 
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this  would  have  produced  an  excessive  amount  of  output  data  and  the  values 
were  therefore  printed  every  MS  or  every  tenth  time  interval.  The  theo¬ 
retical  curves  shown  in  the  graphical  presentation  of  results  were  plotted 
using  points  every  five  milliseconds. 

One  inherent  difficulty  with  the  Runge-Kutta  method,  as  pointed  out 
by  Newmark  and  Chan  (j>J  ,  is  that  "the  amplitude  will  gradually  damp  it¬ 
self  out  even  in  the  case  of  an  undamped  system,  and  the  method  is  there¬ 
fore  not  desirable  for  a  long  period  of  time."  No  explanation  was  given  as 
to  what  constituted  "a  long  period  of  time”;  however,  the  time  limit  imposed 
in  this  problem  was  undoubtedly  short  enough  to  prevent  this  from  affecting 
the  results. 

10.  Perhaps  the  greatest  source  of  error  involved  in  the  experiment  was 
the  reduction  of  the  experimental  data.  The  data  reduction  involved  first 
of  all  a  determination  of  a  zero  reference  level  for  the  various  traces. 

This  was  accomplished  by  using  the  edge  of  the  oscillograph  record  as  a 
reference  and  then  comparing  this  with  the  1000  cycle  trace  for  straight¬ 
ness.  In  all  cases  the  straight  line  formed  by  the  1000  cycle  trace  agreed 
with  the  paper  edge  reference.  Amplitudes  for  the  various  parameters  were 
measured  from  the  respective  reference  lines  using  a  six  inch  steel  scale. 

A  magnifying  glass  was  used  for  greater  accuracy.  The  scale  was  read  to 
the  nearest  0.01". 

As  noted  earlier,  the  sensitivity  of  the  velocity  meters  varied  ±  3% 
during  a  4"  sweep.  The  response  of  the  accelerometers  was  good;  the  only 
drawback  in  their  use  was  the  high  frequency  pick-up  in  the  first  30  MS. 

This  completely  obliterated  the  trace  for  the  first  30  or  40  MS,  prevent¬ 
ing  any  real  comparison  of  the  compatibility  between  velocity  and  accelera¬ 
tion  traces. 
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The  integration  of.  the  acceleration  vs time  plots  was  accomplished 
by  using  a  K  &  E  4236  Compensating  Polar  Planimeter,  serial  number  6098. 

It  is  estimated  that  an  error  of  about  2%  was  possible  using  this  instru¬ 
ment  since  the  areas  Involved  were  of  the  order  0.1  sq.  in. 

Since  the  values  obtained  in  the  data  reduction  were  accurate  to  about 
two  or  three  significant  figures,  the  values  put  into  the  computer  were 
limited  to  the  same  precision,  which  was  also  utilized  Its  presenting  the 
graphical  comparison  of  results. 

11.  It  is  to  be  noted  that  the  parameters  chosen  for  the  comparison  of  the 
results  could  have  been  force  vs.  deflection,  acceleration  vs.  time,  or  some 
other  comparison.  The  velocity  .vs.  time  comparison  was  chosen  because  it 
offered  the  quickest  comparison  with  the  least  manipulation  of  numbers.  The 
other  comparisons  mentioned  above,  would  have  involved  extensive  reduction, 
which  would  have  introduced  the  possibility  of  increased  error. 
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7.  Conclusions  and  Recommendations 


The  gratify  inglyr  close  comparisons  between  e xperimentally  observed 
behavior  and  theoretically  predicted  behavior  of  the  three  mass  structure 
considered  herein  permits  the  conclusion  that  the  response  of  similar 
structures  may  be  predicted  by  mathematical  analysis  using  digital  computer 
calculations  and  involving  force-dc fleet ion  relations  for  the  elastic  and 
plastic  regimes  which  involve  prior  deformation  history  in  the  manner  des¬ 
cribed  in  Appendix  I.  The  analysis  in  the  present  case  was  complicated  by 
a  feature  of  the  experimental  set-up  which  would  probably  have  no  counter¬ 
part  in  actual  practice.  There  should  be  no  essential  difficulty  in  analys¬ 
ing  cases  where  the  excitation  ccn.es  fv on  lateral  forces  having  specified 
variations  with  time  rather  char,  fr u  a  velocity  step  of  one  of  masses, 
although  the  latter  case  is  certainly  of  primary  interest  In  considering 
response  to  explosive  loadings.  Similarly,  there  should  be  no  difficulty 
in  incorporating  further  refinements  regarding  "large"  geometry  in  the 
analysis  should  it  be  desired  to  follow  the  deflection  histoiv  into  ranges 
where  the  present  restrictions  to  "small"  geometry  might  not  be  sufficiently 
accurate. 

The  present  experiments  d*  w  t ^ at  for  the  higher  response  vtl  clties 
that  accompany  the  larger  input  velocity  steps,  the  w. !ue  of  cr  chat  gives 
the  best  theoretical  prediction  of  the  response  also  increases,.  For  reasons 
cited  by  Kurzenhauser  M  ,no  attempt  was  made  in  this  thesis  f;  r  aix  a 
quantitative  evaluation  of  this  rehit  ion. 

The  predicted  response  depends  greatly  on  the  details  of  the.  load- 
deflection  relations  which  depend,  in  turn,  upon  a  knowledge  of  the  dynamic 
behavior  of  the  material  as  it  is  loaded  and  unloaded  in  both  the  elastic 
and  plastic  regime.  It  would  be  of  great  advantage  to  have  better  informa¬ 
tion  than  seems  to  be  available  presently  abaut  structural  materials  loaded 
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in  such  a  manner;  including  the  variability  of  c r  with  strain-rate. 

Future  experimental  investigations  should  employ  apparatus  in  which 
a  closer  control  over  or  knowledge  of  column  fixity  is  available.  Also, 
if  a  ballistic  pendulum  type  of  apparatus  is  to  be  used,  it  is  clear  that 
the  masses  should  be  hanging  freely  before  impact  so  as  to  avoid  the  compli 
cations  introduced  in  the  present  experiments  by  the  presence  of  the  hold¬ 
back  cable.  Finally,  it  would  clearly  be  quite  advantageous  to  have  a 
single  instrumentation  capable  of  reporting  the  response  of  the  structure 
beyond  the  100  MS  cut-off  without  the  necessity  of  augmenting  the  response 
with  a  second  instrumentation  as  was  necessary  in  this  investigation,  where 
in  velocity  meter  data  was  extended  by  use  of  integrated  accelerometer 
data.  - 


59 


Bibliography 

1.  Alfred  Kurzenhauser,  Portal  Frame  Under  Impact  Loading,  Thesis, 
U.  S.  Naval  Postgraduate  School,  Monterey,  California,  1961, 

2.  J.  M.  Massard  and  R.  A.  Collins,  The  Engineering  Behavior  of 
Structural  Metals  Under  Slow  and  Rapid  Loading,  Structural 
Research  Series  No,  161,  University  of  Illinois,  October  1958. 

3.  A*  H.  Keil,  The  Response  of  Ships  to  Underwater  Explosions, 
Structural  Mechanics  Laboratory,  Report  1576,  David  Taylor 
Model  Basin,  November  1961. 

4.  G.  L.  O'Hara,  Shock  Spectra  and  Design  Shock  Spectra, 

Structures  Branch  Mechanics  Division,  U.  S.  Naval  Research 
Laboratory,  November  12,  1959. 

5.  N.  M.  Newmark  and  S.  P.  Chan,  A  Comparison  of  Numerical  Methods 
for  Analyzing  the  Dynamic  Response  of  Structures,  Structural 
Research  Series  No.  36,  University  of  Illinois,  October  1952. 

6.  S.  Timoshenko,  Strength  of  Materials  -  Part  I,  3rd  Edition, 
p.  267,  D.  Von  Nostrand  Co.,  Inc.,  New  York,  N.  Y. ,  1955. 

7 .  Alcoa  Aluminum  Handbook,  Aluminum  Company  of  America, 

1957. 

8.  W.  E.  Milne,  Numerical  Solution  of  Differential  Equations, 

John  Wiley  &  Sons,  1953. 


60 


61 


k>3 


ure 


ACCELE  RO  METER'S 

;i  it 


V5L0CI"  Y 
METERS 

”  ’1  ’• 

.  1 1  .  , 


a? 


I  IOOO 
~1  CYCLE 


CONTROL  'Jl. 


6  V 


RECORDS-  f 
CXILLOGR-APH! 

i 


Figure  4 


Figure  5a 


x3n. 

i  => 


x. 


i  ^ 

X-14. 1 
yJ 


X- 


M 


>Wt 


T 

R 


><, 
•  A 


Figure  5b 


64 


AK'StLT  BFTHZVll  / 
VERTICAL  Ai.’D  I 
SUrPORTINS-  4  | 
WIRE - 

7  I 
/  I 


R 


\ 


Tv  /' 

A  / 


v 


y 

A 


4->lHf  x-  <r 


V 

V 


Figure  6a 


Figure  6b 


ll: 

J 

vy 

mj 

U 

14 

o 

Ll 


K- 


X*  >1 


\  / 

4- 


/-K 


,.X- 


/] 


I  / 
!/ 


,V  L/riV"  '  CEELF  C'UC'W  (IK.) 


A 


Figure  7 


65 


66 


Total  Momentum  (lb-sec) 


56 


j 


26 


24 

0  20  40  60 

Time  (MS) 


80  100 


Figure  9 

Plot  of  total  momentum  vs,  time  for  an  initial  velocity  input 
(Step)  of  1.82  ft/sec.  The  symbol  O  denotes  experimental 
points  and  the  solid  curve  represents  theoretical  values  for 
an  initial  anvil  displacement  of  0.7  ft. 
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AFFENDTX  I 


FLOW  CHARTS  AND  COMPUTER  PROGRAM 
This  section  contains  a  very  general  flow  chart  for  the  main 
program,  a  more  detailed  flow  chart  of  the  force  versus  deflection 
routine  which  was  the  controlling  routine  of  the  program,  and  the  actual 
computer  program.  The  force  versus  relative  deflection  routine  permitted 
the  evaluation  of  the  spring  forces  as  the  mathematical  model  experienced 
elastic  -  plastic  response.  A  list  of  symbols  is  included,  since  several 
symbol  changes  were  necessary  in  conversion  to  the  Fortran  language*  The 
figure  shown  belo\  and  the  force' deflect  ion  flow  chart  shown  on  page  74 
explain  the  force-deflection  subroutine. 


The  subroutine  is  entered  from  the  main  program  shown  by  the  box  labeled 
HP.  Several  parameters  are  picked  up  such  as  the  relative  displacement  of 
the  columns,  denoted  by  S,  the  abscissa  3  of  the  knee-point,  the  elastic 
and  plastic  spring  constant  r»  k  and  k*,  the  value  for  F  where  the  extension 
of  the  line  represented  by  k*  crosses  the  vertical  or  force  axis  (denoted 
by  A),  and  the  coordinates  (P,G)  of  the  present  location  of  the  origin  of 
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the  force-deflection  plot.  Note  that  k*  is  negative  as  derived  on  page  15 
This  is  indicated  in  the  figure  above.  The  symbol  R  is  the  value  of 
displacement  used  throughout  the  routine.  This  value  is  continually 
corrected  for  the  changing  origin.  The  location  of  the  origin  is  deter¬ 
mined  from  the  previous  history  of  the  response.  The  first  time  that  the 
subroutine  is  entered  the  values  of  P  and  G  are  zero,  that  is  the  origin 
is  at  (0,0).  If  the  relative  displacement  S  is  less  than  the  knee-point 
value  B,  then  the  force  will  be  equal  to  the  elastic  spring  constant  k 
times  the  relative  displacement  R,  where  R=S-PorR=S  since  P  = 
zero.  When  the  value  of  R  builds  up  to  exceed  the  knee-point  value  B, 
a  check  is  made  to  see  if  the  column  is  loading  or  unloading.  The  check 
used  is  a  displacement  check  denoted  by  the  symbol  U,  where  U  is  equal 
to  the  present  relative  displacement  S  minus  the  previous  relative  dis¬ 
placement  SS.  As  long  as  the  present  displacement  is  larger,  the  forcing 
function  is  equal  to  A  +  k*R.  The  above  routine  is  followed  until  the 
columns  begin  to  unload.  The  present  value  of  relative  displacement 
becomes  smaller  than  the  previous  value  and  a  new  origin  is  computed 
from  the  relations:  P  =  0  +  H  -  B  and  G  =  k*(P) 

where  0  is  the  horizontal  value  of  the  previous  origin  and  H  is  the  last 
recorded  value  of  R.  The  routine  loops  back  to  R  =  S  -  P  where  R  will 
now  be  less  than  the  knee-point  value  B.  The  value  for  force  will  be 
computed  from  F  =  G  +  kR  or  for  values  greater  than  BbyF=A+G+  k*R. 
The  above  described  routine  continues  until  the  problem  ends  and  a  new 
step  input  is  introduced  whereupon  all  parameters  are  initialized.  The 
same  type  of  routine  is  also  included  for  the  case  where  the  relative  dis¬ 
placement  or  forcing  function  or  both  become  negative.  This  can  be  seen 
by  following  the  same  type  of  analysis  used  above. 
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The  force-deflection  routine  is  found  inside  the  subroutine  labeled 


SUBROUTINE  DERIV  statement  number  77  of  the  program  listing,. 
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List  of  Symbols 


STEP  -  the  initial  input  velocity  to  anvil  (ft/sec) 


V 


time  increment  (sec) 


XE  -  time  limit  of  problem  (sec) 

DE  -  end  value  of  step  input  (ft/sec) 

DD  -  initial  anvil  displacement  [a]*  (ft) 

CHG  -  correction  factor  to  allow  for  variable  initial  anvil 
displacement 

DC  -  ^anvil  cable  length  -  [dd]  a  (ft  ) 

YD.  -  displacement  corresponding  to  ultimate  bending  moment 


STR^  -  plastic  spring  constant  [k*]  (f// in") 

FBR.  -  intersection  of  plastic  force-deflection  line  wit 
vertical  axis  (lbs) 

ANS  -  output  printing  interval  (sec) 

SP(I)-  step  input  values  for  I  different  values  (ft/sec) 


-  independent  variables 

-  dependent  variables 

F  -  forcing  function  in  columns  (lbs) 


S 


present  relative  displacement 


SS  -  previous  relative  displacement 


X 


time  variable 


TAB  -  time  expired  between  print  intervals 


0 


last  value  for  relative  displacement  of  origin  of  force- 
deflection  plot 


H 


last  value  for  force  origin  of  force-deflection  plot 


^Symbols  in  refer  to  previously  defined  terms. 
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G  -  new  force  coordinate  for  force-deflect  Lon  plot 

P  -  new  relative  displacement  coordinate  for  force-deflection 
plat 

R  -  Relative  deflection  of  the  columns 
U  -  S  — SS 

SUM  -  total  momentum  of  the  structure 
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Main  Flow  Chart 


Statement  // 
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**  Statement  numbers  are  sequential  starting  with 
PROGRAM  KENBOB  as  statement  number  one. 
Continued  statements  count  as  one  statement. 
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Force-Deflection  Flow  Chart 
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*  P  =  G  =  Zero,  initially  upon 
first  entry  from  main  program. 


Computer  Program  for  Steel  Columns* 
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*The  program  for  aluminum  columns  is  identical  except  for 
slightly  different  values  for  and  M^. 


75 


G(  2  ) 
H(  1  ) 
H(  2  ) 


35 

STEP  =  SP( I ) 

I  F ( STEP-DE ) 35,  16 

,  16 

I F ( SEN SE  SWITCH 

1 ) 33, 34 

33  PRINT  3  1  , STEP 

310  FORMAT (7H  STEP 

=  F5.2/// 

1G1H 

TIME 

X(  1  ) 

1  V  (  1  ) 

F  (  1  ) 

F  ( 2  ) 

S{  1  ) 

S  (  2  )  ) 

GO  TO  32 

34 

30 

PRINT  30, STEP 

0  FORMAT ( /H  STEP 

=  F5.2// / 

10  1H 

TIME 

SUM 

32 

1  DD 

V(  2) 

V  ( 3  ) 

A  ( 2  ) 

A  (  3  )  ) 

Y ( 4 )  =  STEP 

10  X  =  X  +  V 

CALL  RKUTT  A ( 6 , X , Y , V ) 

D(5)  =  ( F ( 1 )-F(2) 1/0.743 
0(6)  =  F ( 2 ) / 0 . 7  1 7 

SUM  =  18. 90* Y ( 4 )  +  0 . 743 * Y ( 5  )  +  0.717*Y(6) 

00  =  DD  +  CHG 

DC  =  64.00  -  00**2 

TAB  =  TAB  +  V 

I F ( TAR-ANS  )  12,18,18 

18  TAB  =  0. 

I F ( SENSE  SWITCH  1  )13,14 

13  PRINT  19,X,Y(  I  ) ,Y(4) ,F( 1  ) ,F(2) ,S(  1  > ,S(2) 

GO  TO  19 

14  PRINT  19, X, SUM  ,00  ,  Y  ( 5  )  ,Y(6)  ,0(5)  ,0(6) 

19  FORMAT  (  1  PE  15.4  ,  6E  15.3) 

12  I F ( X-XE  )  10,  15,  15 

16  1  =  0 


PAUSE 

. OREAD.  17.,CHG.  ,V,XE,DE  ,00,  YD  (  1  )*  YD  (2)  ,  SL  (  1  )  ,  SL  (2  )  ,  STR(  1  )  ,  STR(2  )  , 

'  ■  1 FBR ( 1 ) ,FBR(2 ) , ANS,SP( 1 ) ,  SP ( 2 ) , SP ( 3 ) , SP ( 4 ) , SP ( 5 ) , SP ( 6 ) 

17  FORMAT ( 10E8.2) 

DC  =  64.00  -  DD**2 

GO  TO  15 

END 

SUBROUTINE  HER  I  V  (D,Y,X) 

GDI  MENS  I  ON  D(6'),Y(6),R(2),S(2),F(2)»0(2),G(2),H(2),P(2),U(2), 
1YD(2)  ,SL(2)  ,  STRt  2)  ,Ff,R(2  )  ,  SS  (  2  ) 

-  -  COMMON  - YO , SL , ST? , FBR , F , S  S , 0  »  G, H, S, DC , DC 

0(1)  =  Y ( 4  ) 

D( 2)  =  Y( 5) 

D  (  3  )  =  Y  (  6  ) 

S (  1  )  =  ( Y{  1)  -  Y ( 2  )  )  *  12. 

S ( 2  )  =  ( Y ( 2 )  -  Y ( 3  )  )  *  12. 

U(  1 )  =  S(  1 )  -  SS (  1  ) 

U  (  2  )  =  S  (  2  >  -  S  S  i  2  1 


9  R (  1  )  =  S(  1  )  -  0(1) 

I F ( R (  1)-Y0( 1  )  )  1,3,3 

1  I F ( R {  1  )  +  YD{  1  )  )5,5,2 

2  FID  =  G(l)  +  SL(l)*Rll) 

SS( 1 )  =  S(l) 

GO  TO  11 

3  IF { U ( 1 ) )  7,4,4 

4  F(  1 )  =  FBR(l)  +  G(  1  )  +  STR ( 1  )  *R  (  1  ) 
H(  1  )  =  R  (  1  ) 

SS (  1 )  =  S (  1  ) 

GO  TO  1 1 

5  IF ( U ( 1 ) )  6,6,8 

6  Ft  1  )  =  -FBR(l)  +  G (  1 )  +  STR { 1 ) *R (  1  ) 
H(  1  )  =  R  (  1  ) 

SS (  I  )  =  S(  I) 

GO  TO  1  1 

l  P(1)  =  0(1)  +  H(l)  -  YD(1) 

G(1 )  =  STR ( I ) *P (  1  ) 

0(1)  = 

GO  TO  9 

8  P(l)  =  0(1)  +  H(l)  +  YD(1) 

G (  1 )  =  STR {  1  )*P(  1  ) 

0(1)  =  P(  1  ) 

GO  TO  9 

l 1  R ( 2 )  =  S ( 2 )  -  0(2) 

IF ( R( 2 )-YD( 2 ) ) 2  1 , 23, 23 

21  IF(R(2)+YD(2) )25»25,22 

22  F  (  2  )  =  G(  2  )  +  SL  (  2  )  *R  (  2  ) 

SS ( 2 )  =  S (  2 ) 

GO  TO  29 

23  IF ( U( 2 ) )  27,24,24 
24  F( 2 )  =  F8R ( 2 )  +  G(2)  +  STR(2)*R(2) 
H(  2  )  =  R(  2  ) 

SS  (  2  )  =  S  (  2 ) 

GO  TO  29 

25  IF ( U ( 2 ) )  26,26,28 

26  F( 2 )  =  -FBR ( 2 )  +  G(2)  +  STR(2)*R(2) 
H(2)  =  R ( 2 ) 

SS  (  2 )  =  S ( 2 ) 

GO  TO  29 

27  P ( 2 )  =  0(2)  +  H( 2 )  -  YD ( 2 ) 

G ( 2 )  =  STR ( 2 ) *P ( 2 ) 

0(2)  =  P ( 2 ) 

GO  TO  11 

28  P ( 2 )  =  0(2)  +  H( 2 )  +  YD ( 2 ) 

G ( 2 )  =  STR ( 2 ) *P ( 2 ) 


t' 


0(2)  =  P  (  2 ) 

GO  TO  11 

29  CC  =  DO  +  Y(  1) 

YX  =  8.00  -  SQRTF(DC-CC**2) 

CK  =  8.00  -  YX 

ODUT  =C^-F (  1  )-(CC/CK)*656.*(  l.+( Y(4)**2+YV**2)/(32.20*CK  )  )  )/(  18.9 

1 +(656./ 32.20 )* (CC/CK ) **2 ) 

0(5)  =  ( F (  1  )  -F ( 2  )  J/0.743 

0(6)  =  F ( 2  > / 0 . 7  1 7  . 

ENO 

SUBROUT  TNi  >  RKUTT  A  (NUMBER  ,  XVAR, YVARS , STEP ) 

DIMENSION  YVARS (30 ) » AK(4  »30) ,  DY ( 30 ) , YC { 30 ) , C ( 3 ) 

C(  1  >=0.5 
C( 2)=0.5 
C  (  3  )  =  1 . 0 

CALL  OER I  V ( OY ,  YVARS, XVAR) 

DO  1  J=1, NUMBER 

1  AK (  1, J  )=STEP*DY( J) 

DO  2  1=2,4 

XC  =  XV  AR  +C ( I-  1  )*STEP 
DO  3  J=l, NUMBER 

3  YC ( J ) =  YVARS(J)  +  C ( 1-  1 )  * AK ( I- 1 , J ) 

CALL  DERIV(DY,YC»XC) 

DO  2  J=l, NUMBER 

2  AK ( I , J  )  =STEP  *0Y{  J  ) 

DO  4  J=l, NUMBER 

4  YVARS ( J )  =  YVARS( J )  +  ( AK( 1 , J ) +2 . * AK ( 2 , J ) +2 . * AK ( 3 , J ) +  AK ( 4 , J )  ) /6 . 0 
RETURN 

END 

END 


APPENDIX  II 


REDUCTION  OF  DATA 

This  section  contains  a  complete  set  of  calculations  for  one  run. 
Also  included  is  a  tracing  of  the  oscillograph  record  for  the  same  run. 

A  photograph  of  the  actual  record  was  not  possible  due  to  the  lack  of 
contrast  between  the  trace  and  the  paper.  It  is  to  be  noted  that  the 
1000  cycle  trace,  mentioned  earlier,  is  not  shown.  This  particular 
trace  was  omitted  to  prevent  cluttering  of  the  other  tracings  and  also 
because  of  the  tedious  tracing  job  required  to  accurately  reproduce  this 
trace. 

The  data  reduction  consisted  of  the  following  steps: 

1.  Establishing  reference  lines  on  oscillograph  record 

2.  Measuring  trace  amplitudes  for  velocity  parameters 

3.  Integrating  the  acceleration  vs.  time  trace 

4.  Computing  conversion  factors 

5.  Converting  measured  parameters  to  proper  units. 

The  run  in  question  was  for  steel  columns  and  a  step  input  of  5.02 
ft/sec. 

Computation  of  the  necessary  spring  constants,  and  other  input  para¬ 
meters  are  also  shown. 
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1 

Trace  of  Oscillograph  Record 


The  preceding  figure  shows  traces  of  the  oscillograph  record. 
The  symbols  V  denote  velocity  and  A  denote  acceleration.  Direct 
measurements  from  the  figure  gives  the  following  values  in  inchess 


Time 

V(1) 

V  (2) 

V(3) 

0 

0.88 

0 

0 

10 

* 

.18 

0 

20 

* 

.59 

.08 

30 

* 

.75 

.37 

40 

* 

.75 

.85 

50 

* 

.86 

1.24 

60 

* 

1.14 

1.40 

70 

* 

1.39 

1.37 

^values  not  measured 
Calibration  pulse  measurement s s 

V(l)  -  0.720  inches 

% 

V(2)  -  0.785  inches 
V(3)  -  0.820  inches 
A(2)  -  1.24  inches 
A(3)  -  1.12  inches 

Velocity  meter  sensitivities  and  accelerometer  calibration  constants 

V(l)  -  146  mv/fps 
V(2)  -  167  mv/fps 
V (3)  -  151  mv/fps 
**A(2)  -  70.7  g 
**A ( 3 )  -  78.0  g 


**  values  obtained  by  the  David  Taylor  Model  Basin  using  calibration 
resistance. 
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Instrument  settings : 

a.  Calibration 

Amplifier  gain  set  at  30  for  accelerometers,  paper  speed 
r,et  at  4"/sec 

Battery  voltage  for  velocity  meters  set  at  0,6  volts 

b .  Run 

Amplifier  gain  set  at  5  for  accelerometers,  paper  speed 
set  at  40M/sec 

Combination  of  the  above  information  gives  the  following  conversion 
factors : 


a.  for  velocity  conversion 


VpaCTO  R 


CAL  I  6  .  YQLTAQ  £  C  v) 
SENSITIVITY  (y/tys) 


PULSE  MX 


Sample; 


V 


I  FACTOR  — 


0-6  V  #  _ !__ 

0J46  *  0-720  10 


-  517/ 


VXrtCT 5  7 

v3  tact.e  —  4.  84- 


b.  for  acceleration  factors 


A  "AC  TO  K 


CALlbS-  CQMS~/^JT 
CALll/j,  PULS).  HT. 


Run  GvAiim 

CALl6„&/U  IV) 


Sample 


A , 


r  ACTOR, 


 (70.  7)  (52. 2)  ^-Ca 


.24- 


n.. 


30 


=  506 


R/\C  TOC 


3  74- 


c.  conversion  for  integration  constant 


vc- 

Vat 


FflCTotf 

P  ACT  O  (2. 


Vac- 

V^C-  PrtCTOR 


k  ac  *r  o  *  T  IM  t  t  actOR.  *  A  h?  E  A 
306(%  ^  025  ■  AREA(/) 

7  •  Ak'r  A  C  'V'rc  j 

9.  5  5  •  AREA  (  fVW  ) 
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Areas  under  acceleration  trace: 


Time  range  (MS) 

Area  V(2)  (in^) 

Area  V(3) 

70-80 

0.01 

-0.04 

70-90 

-0.14 

-0.06 

70-100 

-0.32 

-0.14 

From  the  above  information: 


Time  (MS) 

V(l)  (ft/sec) 

V(2)  (ft/ sec) 

V(3)  (ft/ 

0 

5.02 

0 

0 

10 

k 

0.82 

0 

20 

* 

2.70 

0.39 

30 

* 

3.43 

1.79 

40 

k 

3.43 

4.12 

50 

* 

3.92 

6.00 

60 

*  5 

5.21 

6.78 

70 

* 

6.35 

6.63 

80 

k 

6.43 

6.26 

90 

* 

5.28 

6.07 

100 

k 

3.90 

5.32 

*Not  measured 


Evaluation  of  spring  constants  -  steel 

Column  thickness  *  h  ^  0.128’1 

Column  width  =  b  =  1.95" 

Effective  length  =  1  s  90." 

ef  f 

g 

■^modulus  of  elasticity  =  E  ~  29.4  x  10  psi 

Weight  on  upper  columns  ^  =r  23.1  lbs 

Weight  on  lower  columns  =  Q  =  47.0  lbs 

*Value  obtained  from  actual  tension  test  of  sample 
column.  See  Appendix  III 
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I  - 
Q*  ~ 

K,  = 


b  h' 
I  2" 


_  (/  SS)  (  -  I  i£  X 


i  2 


_27T  jEI  _  2  TT^luO  4-)(  54  I  ) 


-  4’ 

34/  x  »  o  /  n . 


lbs. 


8  I 


24  t  X 


0  3 


1  - 

Q,  I 

iu  4  )  (  29  4-X  34  I) 

1  4  /  0  n 

9 

Q*J 

9  3 

2.44  Oj 

K, 

K2 


-  32  4-  '  D/i* 

=.32  1  •*■/>» 


K*  = 


Qj 

i 


41  0 


5.  2  2 


'V- 


I  r\  , 


/</  = 


2.51  'Vi 


/h 


Evaluation  of  force-deflection  knee  points 


(a-.)  er  =.  <40,  000 

cr  b  hz  -  C4-y,ooa)(l.9sX°  -  I27&ir\-lb 


«,  X,^(  - 


324  Xu,  - 


4  Mu  "  Qi  X41 
_ 

ai8  —  47.0  Xu., 


X  U2. 

A 


1 2.  7© 

A9  <b  3 

—  O.  4  3  1 

1  O 

1273 

29  6  & 

-  0-420 

/  fo 

4  Mu 

1  i  7  B 

Jl 

9 

142.  It 
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(b) 

cr  ~ 

So,  ooo 

A*  «- 

,  — 

o,  50 

1  *  1 

X  - 

O.  5  Sfe 

/  n 

A  = 

/  7  7  /  b 

CO 

CT  rr 

Go,  ooo 

X  u  i  — 

O  fe-4  6 

//n 

X  LA,  Z  “ 

O.  64  5 

/  >i . 

A  = 

2/3  1  L, 

Using 

the  same 

procedure 

for  3003-H14 

a luminum: 

K,  n 

><Z  - 

k;  = 
k;  = 

103  1 b /m 

* 

112  ,b/  i* 

-  5.2  |  lb/,n, 

-  2  53 

WH  fe  k  F-  •  h  - 
b  - 

<*.  - 

Km  e  E. 

Po  1  M)  T  £ 

• 

CO 

cr  - 

2  0;  OOO 

pS  c 

Xu,  - 

0.  fe  19 

/  1A 

Xu*. 

0  6  1  (6 

1  lf\ 

A  - 

70  ft 

/  b  . 

(k) 

cr-  - 

50,  OOO 

r* 

Xu,  - 

0.^28 

Am-,  . 

0.42S 

/  ^ 

A  - 

1  O  6 

i  b. 

O,  US 

l .  9  S " 

=  ft  15  " 

G 

/  o  o  y  i  o  pi 
22  , 


t  t 


lvalue  taken  from  Alcoa  Aluminum  Handbook  for  3003-H14  Aluminum 
reference  £?j  . 
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APPENDIX  III 


COLUMN  MATERIALS,  TENSION  TESTS,  AND  EFFECTIVE  MASSES 


Column  Materials 

As  noted  earlier,  two  different  materials  were  utilized  in  the  test 
columns;  hot  rolled  steel,  with  properties  similar  to  those  of  the  steel 
used  by  Kurzenhauser,  and  3003  -  H14  aluminum.  The  aluminum  sheet  used 
had  a  yield  strength  of  21,000  psi  and  an  ultimate  strength  of  22,000 
psi,  as  listed  by  the  Alcoa  Aluminum  Handbook  7  .  A  6061  T6  aluminum 

used  by  Massard  and  Collins  2  ,  would  have  been  more  representative 

of  the  aluminum  used  in  structural  work;  however,  this  was  not  available 
at  the  time  the  experimentation  was  conducted. 

a 

Tension  Tests 

Three  samples  of  each  of  the  column  materials  used  were  selected  and 
shaped  into  standard  tensile  test  specimens  with  a  2"  gage  length  to  meet 
ASTM  Standard  E8-54T.  A  Baldwin  Universal  Testing  machine,  along  with  a 
Riehle  extensometer  were  utilized  in  making  the  tension  tests.  The  results 
are  tabulated  in  Table  1  below: 


Table  1 


Steel 

g 

Modulus  of  Elasticity  29.4  x  10  psi 


Aluminum 


Yield  Point 


40,800  psi 


20,300  psi 


Ultimate  Strength 


59,300  psi 


22,250  psi 
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Effective  Masses 


A  very  simple  approximation  of  the  effective  weights  of  the  structure 
was  made  by  including  one  half  of  the  weight  of  two  supporting  columns 
with  the  upper  mass  and  the  weight  of  two  full  supporting  columns  with 
the  center  mass.  The  remaining  two  halves  of  the  lower  columns  were  in¬ 
cluded  with  the  anvil  mass.  The  figure  below  shows  this  procedure. 


M,  n-M-CTlVi: 


i  -  - 


lH'liCTIVE 


-X- 


r  I  I  I  C  T  i  V  E 


I 


M 


V 


The  values  of  the  above  effective  masses  are  listed  in  table  2 


2 

The  units  are  lb-sec  /ft 


Table  2 


Mass  Number 


Steel  effective 


A1  effective 


1 


18.9 


18.9 


2 


0.743 


0.711 


3 


0.717 


0.705 
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APPENDIX  IV 


FURTHER  PROGRAM  MODIFICATIONS  CONSIDERED 

Two  further  modifications  to  the  computer  solution  were  considered. 
One  modification  was  the  determination  of  the  effect  of  vertical  accelera¬ 
tion  upon  the  structured  columns*  The  fact  that  the  anvil  experiences 
a  vertical  acceleration  upon  impact,  produces  a  force  in  addition  to  the 

structure's  static  weight*  This  vertical  acceleration  was  obtained  for 

2 

several  runs  and  found  not  to  exceed  about  5ft/sec  .  This  would  have 
resulted  in  a  variable  correction  to  the  effective  value  of  Q  of  not 
more  than  167o,  and  the  Q/Qe  correction  was  itself  rather  small,  of  the 

order  of  2%.  Thus  it  was  not  consideted  of  sufficient  importance  to 

* 

make  this  correction  which  would  have  greatly  complicated  the  theoretical 
analysis  and  would  have  resulted  in  considerably  more  cumbersome  dif¬ 
ferential  equations  to  work  with. 

The  second  modification  was  an  attempt  to  force  the  theoretical  and 
experimental  total  momentum  vs.  time  curves  to  agree  exactly,  by  varying 
the  value  of  initial  anvil  displacement  during  the  run  in  such  a  manner 
as  to  retain  the  proper  momentum  relation.  This  was  tried  for  several 
runs,  resulting  in  perfect  agreement  in  the  plots  of  total  momentum  vs. 
time;  however,  practically  no  change  in  the  velocity  response  was  noted. 
This  modification  was  likewise  not  included  since  the  constant  value  of 
a  =  0.7'  produced  essentially  the  same  results. 


88 


APPENDIX  V 


TYPICAL  COMPUTER  OUTPUT 

The  tabulation  of  theoretical  values  for  a  typical  run  shows  the 
time  in  MS,  velocities  V(2)  and  V(3)  in  ft/sec,  and  the  total  momentum 
in  lb- sec  designated  by  the  term  SUM. 
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Computer  Solution  STEP  =  5.02 
(T -  40,000  psi  Steel 


TIME 

1 .  0000E-03 
2 . 0000E-03 
3. 0000E-03 
4 .0000E-03 
5.0000E-03 
6.0000E-03 
7.0000E-03 
8 .  0000F-03 
9.0000E-03 
1 . 0000E-02 
1  .  1  000E-02 
1 .2000E-02 
1 .3000E-02 
I  .4000E-02 
1 .5000E-02 
1 .6000L-02 
1  .  7000E-02 
1 . 8000E-02 
1 . 900 OF- 02 
2 . 0000E-02 

2.  1  000E-02 
2 . 2000E-02 
2. 3000E-02 
2.4000E-02 
2.  .5000E-02 
2.6000E-02 
2 . 7000E- 02 
2.8000E-02 
2.9000E-02 
3 . 0000E-02 
3. 1 000E-02 
3 . 2000E-02 
3 . 3000E-02 
3.4000E-02 
3. 5000E-02 
3. 6000E-02 
3.7000E-02 
3.8000E-0? 
3. 9000 E- 02 
4 . 0000E- 02 
4 . 1 000E-02 
4 . 2000E-02 
4. 3000 E- 02 
4.4000E-02 
4 . 5000E-02 
4 . 6000E- 02 
4.7000E-02 
4 .8000E-02 
4 . 9000E-02 
5. 0000 E- 02 
5. 1 000E-02 
b .2000E-02 
5 . 3000E-02 
5.4000E-02 
5.5000E-02 
5.6000F-0? 
5. 7000 E- 02 
5 . 8000E-02 
5. 9000E-02 
6 . 0000E-02 


SUM 

9 . 4  8  1  E  +  0 1 

9.475L-  +  01 
9.469E+01 
V.462L+01 
9.456E  +  01 
9.449E+01 

9 . 4  4  3  E  +  0 1 
9.437E+01 

9.4  30E  +  01 
V.423E+01 

9.4  17E  +  01 

9.4  10E  +  0I 
9.408E+01 

9. 3  y  7  E h  0 1 

9 . 3  9  0  E  +  0 1 
9. 3 83 E +01 
9 . 3  7  6 1  +  0 1 
9. 3>6  9  If  +  0  I 

9 . 3  6  2  F  +,  0 1 
9.  3bbE  +  01 
9.348E+01 
9.34  IE  +  01 
9. 334 E+ 01 
9. 327 E +01 
9. 320E+01 
9.313E+01 
9.306E+01 
9.2981 +01 
9.291 6+01 

9 . 2  84t  +  0 1 

9. 2  76F  +  01 
9.269i_ +  01 

9 . 2  6  1  E  +  0 1 
9.254E+01 
9.246E+01 
9.239c+01 
9.231E  +  01 
9.223E+0I 
9.2  I  6  r  +  0 1 
9 . 208t +0 1 
9.200E+0) 
9. 193E+01 
V. 1 85  b  +  0 1 
7.  17/E  +  01 
V.  1  69E  +  01 
9.  1  6  1  i:  +  0  1 
9.  1S3E  +  01 
9. I 45t+0 1 
9. 137E+0I 
/.  1  2  9  E  +  0 1 
9.  1  2  1E  +  01 
9. 1  I  3E  + 01 
9.  1  ObF+Ol 
9.096E+0I 
9. 0881: +  01 
9.080F+Q1 
9. 072E+01 
9.063E+OI 
9 . 055E  +  Oil 
9.047F+01 

li 

0 


V  ( 2  ) 

1 .3  12E-02 
S.233E-02 
1  .  I  /2E-0 1 
2.090 E - 0  1 
3.20  7F-0  1 
4 .  b  7  3  E  -  0  1 
6.  lbOE-O  1 
7 . 902E-0 I 
9.645E-0 1 
1  .  134 E+00 
1 .297E  +  00 
1 .45  3E  +  00 
1 . 60  2  E  +  00 
1 . 74 3 t +00 
1 .876 E+00 
2 . 000 E  +  00 
2. 1 14E+00 
2.220E+00 
2.3 16 E+00 
2.402E+00 
2.479E+00 
2.'.»4  7r-  +  00 
2. 606E+00 

2. 6  b  7 1: +  00 
2 .699 E+ 00 
2. 734 E+00 
2 . 763E+00 
2. 985E +00 
2.80 1 E+00 

2.8  1  31: +  00 
2. 822 E  +  00 
2 . 827  E  +  00 
2.831 E+00 
2.8341:1  00 
2.637  F  +  00 
2.84  1E  +  00 
2.84  71:  ♦  00 
2. 656 E+00 
2. 869 E  +  00 
2.88/ E+00 

2.9  10E  +00 
2.939b +00 
2. V  78  E  +  00 
3.0 1 9 E  t 00 
3.0  7 1 E  + 00 
3.131  F ♦ 00 
3 . 200 1:  +  00 
3.278F  » 00 
3. 366 E  +  00 
3. 4  621.  *  00 
3.5681 +00 
3.683E+00 

3 . 6  0  7  F  t  0  0 
3.940E+ 00 
4.080 11  +  00 
4.228F +00 
4.  38  3  f-*  00 
4.54  bl. +00 
4.71  11  +00 
4 .8801+00 


V(  3  ) 

b 

.  98  61-0  6 

9 

.b6  IF— Ob 

4 

.82 / E -04 

1 

.b  I9L-0  3 

3 

.69  IF  -0  3 

7 

.  60m— 0  3 

1 

. 3991-02 

2 

.  3681—02 

3 

.  /bbF-08 

S 

.649F-02 

H 

.  1  3 3E-02 

1 

.  1  28E-0  I 

1 

. b  I6E-0 1 

1 

.98  31-0  1 

2 

.  b  3  b  1-0  1 

3 

.  1  76F-0  1 

3 

.  90HE-0 1 

4 

.  7  851—0  1 

b 

.  6b  7  F— 0  1 

6 

.6  75F-0  1 

7 

.  7881—0  1 

8 

.99  31—0  1 

1 

.02  91.  *  00 

1 

.  16/1.  100 

1 

.  3  1  31-  i  00 

1 

.  4  6  7  E +00 

1 

. 62 8 f  +00 

1 

.  7VbE  +00 

1 

.9681  i  00 

2. 

.  1 4  b E i 00 

2 

.  3 2 b I-  +  00 

2. 

. bO 9  E  + 00 

2, 

.694 F  ♦  0() 

2. 

.  88011  +1)0 

3. 

, 06 6 1.  •  Of) 

3. 

.  2 b 0 E +00 

3. 

.  4  3  21.  <  0  0 

3, 

.61  1  E *  00 

.  7861  «  00 

3 « 

.9b 61  i  00 

4. 

.  1201-t  00 

4  . 

277 li 00 

4-. 

4  2  71:+  0  0 

4. 

5 6 9 1:  i  Of) 

4'. 

. 70 3 E i 00 

4  . 

82 /E+00 

4. 

9  )|  8  F  +  0  0 

b. 

04  9I.-+  00 

5. 

14  51:  +  00 

b  . 

2  3  21-  +  00 

b. 

30  81- +  00 

bl. 

3/61;+  00 

5. 

4  3 1)  E  «  00 

b. 

4  8  3(1  i  00 

V. 

52  4 F t  00 

5. 

b  5  6 1  • 00 

b  . 

b82l.  t  00 

b . 

60  1 1  t  ()() 

b. 

6  141+  00 

5. 

622E+  00 

90 


TIME 

6.  10006-02 
6-20006-02 
6 .  30C0F-C? 

6 . 4  0  0  0  f:  -  0  2 
6.50006-02 
6.6000E-02 
6. 70006-02 
6.80006-02 
6.9000E-02 
7 . 00006-02 
7. 10006-02 
7.20006-02 
7  .30006-02 
7.4000E-02 


7.5000E-0? 
7. 6 000 E- 02 
7 . 7000E-02 
7 . 80006-02 

7.90006- 02 
8. 0000 E- 02 
8-  1 0  0  0  E - 0  2 
8-2000 E- 02 
8 . 3000F-02 
8.4000E-02 
U.5000E-02 
8 . 6  000E-02 
H. 700 Oh- 02 
8.8000E-02 

8- 90006-02 
9.0000E-02 
9. I 000E-02 
9 .20006-02 
9 . 3000E-02 

9- 4000E-0? 
9.50006-02 
9 . 6000E-02 
9. 7000E-02 
9  -  8  0  0  0  E  -  0  2 

9 .90006- 02 
I  .00006-0  1 


SUM 

9.0386 ♦ 01 

;.O  506  4  oi 
7 -02  16  +  0 1 
) .  0  1  5L  +  0! 
9.004E 601 
8. 9936  6 01 
8.9876+01 
8.9786+01 
8.959E  601 
8.9606  +  01- 
8.95  IE  +  01 
8 .9438  +  01 
8.9346+01 
8.9256+01 


8.9  186  +  01 
'■.9066  +  01 
8 . 8  9  7 1  +  0 1 
<>, .  888.6  +  01 
1.879E  +  01 
8 .  8696  +  01 
8.8606+01 
8. 38 lt+01 

0.34 il +01 

6. .  8  32 1-  +  0 1 
13.8226  +  01 
0.3  12L  +  01 
6l.80.3L  +  01 
9. 7936 +01 
8.7836+01 
8.7  75h  +  01 
{..7646  +  01 
. 7S4r+01 
8 . 7  4  4  6  +  U 1 
6 . 7346  601 

8.  7246  +  rU 

•<.7  146  +  01 
d. 7031+ 01 
3.69 3 L +01 
.  6  8  3 1.  +  0 1 
.  6  73 l-  +  01 


V  (2) 

8 . u 5 0 E t  00 
5. 2  2  0 1"  +  0  o 
5. 50 86+00 

5.8526+00 
8.7096+00 
5. 659 E+ 00 
5.9796+00 
6.  129E  +  00 
6.266E+00 
6.350E+00 
6. 4396  +  00 
6.5 1 3E+00 
6. 570E+00 
6.6  1  1E  +  00 


6.634E+00 
6.640E+00 
6.o29E+00 
6.60 1 E+00 
6. 5566  +  00 
6 .4966  +  00 
6.42  16  +  00 
6.  3  >26  +  00 
6 . 2  3  1  F  +  0  0 
6. 1 19 E+00 
5.  79  76  +  00 
8.9666+00 
8.  7  2  9  6  +  00 
5 . 5686+00 
5 .4436+00 
5.2976+00 
5. 1526+00 
5.0096+00 
4 . 8706+00 
4.  7  366  +  00 
u .6096+00 
4.490E+00 
4.33  1F  +  00 
4 . 2826+00 
4 . 1 946+00 
4. 1 196+00 


v(3) 

5.6276+00 
‘>.62  8  6+00 
5.62 76+00 
5.6246+00 
5.62  16  +  00 
5.6  196  +  00 
5.6  186  +  00 
5.6196+00 
5.6236+00 
5 . 6  5GF  +  00 
5.64  16  +  00 
8.6566+00 
5.6776+00 
5. 7026+00 


5.7326+00 
5.7676+00 
5.8066+00 
8.85 16+00 
5. 8996+00 
5.95  1  E  +  CO 
6.0066+00 
6.07>  36  +  00 
6.12  16  +  00 

6.  1306  +  00 
6.2396+00 
6.2  '  76  +  00 
6.5826+00 

7,  .  40  4  E  +  00 
6.45  16  +  00 
6.4936+00 
6.5286+00 
6.5566+00 
6.5756+00 
6.5356+00 

6 . 5  3  5  6  +  0  0 
6.5746+00 
0.55  16  +  00 

6.5  1  76  +  00 
6.4  7  1  6  +  00 
6.4126+00 
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Computer  Solution  STEP  =  5.02 


TIME 

CT  =  50  >000  psi 

SUM 

Steel 

V  { 2  ) 

v  (  ;  > 

3 .9*161  -06 
7 . 36  II  -0'i 

I  . 0000E-03 

9 . 4  8  1 1  +  0 1 

1 .3  12F-02 

2 .0000 E- 03 

9.4731 +01 

5.2  334-02 

3 . 0000E-03 

9.469F+01 

1  .  1  724-0  1 

4 . 87  71  -Oil 
1 .3  194-0  3 

4. 0000 F- 03 

9.462E+01 

2 . 0  70 E-0  1 

5. 0000 E- 03 

9.4  56F  +  0I  • 

3. 20 7 E-0  l 

3.69  11-0  3 

6. 0000 E- 03 

9. 4491: +  01 

4 . 5  73 E-0  1 

7.6081-0  3 

7 . 0000E-03 

/.  4434+01 

6.  I30E-0  1 

1 . 39.91  -  07 

8. 0000E-03 

9. 4374  +  01 

7. 924 E-0  1 

2.3681  -07 

9 . 0000E- 03 

9. 4 30b +01 

9.8764-0  1 

3  .  70  7 1  -07 

1 . 000OE-02 

9.4234+01 

1 . 198E+00 

5.6674-07 

1 . 1 000E-02 

7.4  17E  +  01 

1 . 408F+00 

8.20  1 1 -07 

1 .2000E-02 

9.4  lOE  +  Ol 

1.611 E+00 

1  .  14  7.1  -0 

1 .3000E-0? 

9.404F+01 

1 . 805 E+00 

1  .534) -0 

1 .4000E-02 

9.397L+0I 

1.991 E+00 

2.0321  -0 

1 . 5000E-02 

9. 3901+01 

2. 16  74  +  00 

2.6*t  71-0 

1 . 6000E-02 

'■>.  383E  +  01 

2.333E+0O 

3.  3»i  71-01 

1 .7000E-0? 

9. 377b +01 

2. 4 88 E+00 

*1  .  136*  -0 

1 . 8000E-02 

7. 5 70 E +01 

2.63 1 F+00 

3 . 0  78)  -0 

1 .9000E-02 

-*.3631  +  01 

2. 77)4  F  +  00 

6.1161  -0 

2.0000E-02 

9.356F+01 

2.8854*00 

7.2  721  -0 

2. 1 OOOE-02 

9. 3494  +  01 

2. 994 E  +  00 

8 . 3*i 61  -0 

2 . 2000E-02 

9.^424+01 

3. 092 E*  00 

9.9  3  7)-  — 0  1 

2.3000E-02 

9 . 3  3  5  E  +  0 1 

3.  1  79E  +  00 

1  .  1  *t *)4*  00 

2 . 4000E-02 

9.328.-  +  01 

3.2546*00 

1  .  30  64  )  00 

2„ 5000E- 02 

9.321E+01 

3.32DF+00 

1  .  *)  7)11  )  00 

2.6000E-02 

9.3  l*i E +01 

3. 373 E+00 

1.66 14 ♦ on 

2 . 7000E-02 

7.  3 06 E +01 

3. 4 22 E+00 

1 . 832)  )  00 

2 . 8000E-02 

9 . 299F  +  01 

3.4  60E  *00 

2.0  3  31  *00 

2 . 9000E-02 

9. 29? r +  01 

3. 4 90 E+00 

2.26  11  10(1 

3. O000E-02 

9. 285E  +  01 

3.5144*00 

2 . 4  7  64 ) 00 

3. 1 000E-02 

9.2  77)1  +  01 

3.3  33  4  +  00 

2 .6961  ) 00 

3 . 2000F-02 

V  .  2  7  0 1:  +  0 1 

3.546  4*  00 

2.971  4 • 00 

3.3000E-02 

9. 2624 +01 

3.5374+00 

3. 130) i 00 

3.4000E-02 

V.  2  53 E  +  01 

3.3654* 00 

3.  3801  i  00 

3. 5 OOOE-02 

9.248F+01 

3.3  721  *  00 

3.6  |  ?l  )  00 

3. 6000F-02 

9.240E+01 

3.3794*00 

3.8'l  34  »  0(1 

3. 7000E-02 

9.233L+C1 

3.588  4  »  00 

4 .0734 +00 

3 . 8000E-02 

9. 225 E +01 

3.3994+00 

*)  .30  IE +  00 

3. 9000E-02 

9.2 171+01 

3.6 144+00 

4.5244+00 

4. 0000 E- 02 

9.2104+01 

3. 6 33 E+00 

4.74 34+00 

4  .  1  000E-02 

9.202E+C1 

3. 658E+00 

4.935E+00 

4 . 2000E-02 

9. 1944+01 

3.690E*  00 

5.  16  1  4  +  00 

4 . 3000E-02 

9. 1874+01 

3. 7  30 E  +  00 

5.3384+00 

4 . 4000E-02 

9. 1796+01 

3. 778 E+00 

5.547E+00 

4.5 000E-02 

9 . 17 16+01 

3.835 1;  +  00 

3. 72  8E  +  00 

4. 6 000E-02 

9. 163E+01 

3.903E+00 

3. 894 E*  00 

4. 7000E-02 

?. 1554+01 

3.91(0  4  +  00 

6.03  14 ) 00 

4 . 8000E-02 

9. 147E+01 

4 .069  4  +00 

6.  1  VUE  *  00 

4.9000E-02 

9. 1396+01 

4.  16.94+00 

6.3  324  )  00 

5.0000E-02 

9.1311-  +,0 1 

4.2804+00 

6 . 4 33E  *  00 
6.3664)00 

5.  1  000E-02 

9.  1  231:  +  0  1 

4.4034+00 

5.2000E-02 

9. 1  156  +  01 

4.5374  +00 

6.665E  *  00 

5 . 3000E-02 

9. 1076+01 

4.6  8'2  4  +  00 

6.  732E*  00 

5. 4000E-02 

9. 0994+01 

4.8  3,At+0() 

6 . 82 8 F  )  (10 

5 . 5000E-02 

9.091E+01 

4.9984+00 

6.8934+00 

5 . 6000E-02 

9.083E+01 

5.  1674  +  00 

6. 9*1  84  )  00 

5 . 7000E-02 

9 . 0  7  5 1  +01 

5.34'lt:  t  00 

6.99  34*  00 

5. 8 000E-02 

9.0664+01 

5.5  i;8i;+  oo 

7 . 0  7  91-  (  00 

5 . 9000E-02 

9.0584+01 

5.6974+00 

7.03  74  *00 

6.0000E-02 

9.050E  +  0.1 

5.8  7,6  4+  00 

7 .0  784  *  00 

I  « 

I 


.  J 
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TIKE 

SUK 

V(2) 

V(3) 

6.  J  000F-G2 

✓  .  04  16  +  01 

6 . 5  3  L  +  0  0 

7.0  ME  +  00 

6.2000E-02 

/.  0-586 ■»  01 

6 . 72f)  F*  60 

7.0996+00 

6 . 3000L-02 

t  .  02l>  t  +  0 1 

6.  >9  36  >00 

7  .  10  36  +  00 

6.40006-02 

V  .  0  1  o  F  +  0  1 

6 . 8  5  3  6  +  0  0 

7.  1026  +  UC 

6.50006-02 

V  .007:-  +  01 

6 . 7046+00 

7 . 0  >96  +  00 

6.60006-02 

8. 7  9:>L  +  01 

6.84a  F  i-OO 

7.0  >36  +  00 

6. 70006-02 

.  ;9C6  +  01 

6.  >  7  1  I-  +  C0 

7 .0866  +  0  0 

6.8000F-02 

- .  9  8  2  c  -f  0  1  • 

7.,)8b6i  00 

7.0786+00 

6.9000F-02 

-..9751+01 

7  .  1846  +  00 

7. 0706+00 

7.0000F-02 

8 .  9646  -t  01 

7.26o6+00 

7 .06 .36  +  00 

7. 1 0006-02 

c  .  9  5  5  L  +  0 1 

7.33 16+00 

7.0576+00 

7. 2000 F- 02 

8.9466+01 

7.3796+00 

7.0836+00 

7. 30  001: -  02 

8  .  V  3  7 1  +  0 1 

7 .408  6  +  0  0 

7.050E+00 

7 . 4  0  0  0  E  -  0  2 

8. 9281:  +  01 

7.4 186+00 

7.0496+00 

7.5000F-02 

8 . 9  1 9  c  +  0 1 

7.4096+00 

7.0006+00 

7,6000 F -02 

8.9106+01 

7.3816+00 

7. 0536  +  00 

7.7000t-02 

8.901c+01 

7.3386+00 

7.0586+00 

7. 80006-02 

-.6911+01 

7.2706+00 

7.068F+90 

7.90006-02 

o. 8826+01 

7.  1896+00 

7.0726+00 

8.00006-02 

O.873F+01 

7.0908+00 

7.0816+00 

8.  1  OOOF-G2 

3.8636+01 

6. 9766  +  00 

7. 08  96  +  00 

8.20006-02 

8.8546+01 

6.6486+00 

7.097F+00 

8. 30006-02 

6 . 8  4  4 1.-  +  0 1 

6. 7086+00 

7.  1036  +  00 

8. 40006-02 

8.8356+01 

6.5566+00 

7.  1076  +  00 

8 . 50COE-O? 

6.825++01 

6.394E+00 

7.  108F  +  00 

8. .60006-02 

+••8156  +  01 

6.2246+00 

'7.  1056  +  00 

8.70006-02 

.80-)':  +  0  1 

6 .64  86  +  00 

7.0 >86+00 

8. 8000E-02 

6 . 7966+01 

5.8o86+00 

7.0846+00 

8.9000E-02 

8. 786+ +01 

8. 6856+00 

7.0846+00 

9. 000 0E- 02 

8 . 7  7  6  f:  +  0 1 

5.8016+00 

7.0376+00 

9.  1  000E-U2 

6.7666+ 0 1 

5.3196+00 

7.00  16  +  00 

9. 200 OF- 02 

8.7566+01 

5  .  1 .3 96+  00 

6.9556+00 

9. 3000F-02 

8. 746- +01 

4  .9646  +  00 

6.9006+00 

9.40006-02 

8.7366+01 

4 .795 E+ 00 

6.8546+00 

9. 80006-02 

8.72  5  t-  +  j  1 

4 . 6  3  3  6  +  u  0 

6. 7586+00 

9.6000E-02 

8.7156+01 

4 .4806  +  00 

6.6  ‘:>  96  +  0  0 

9.7000E-U2 

‘<•70  5  »■•  +  0  i 

4 . 3  3  8  h  +  0  0 

6.5696+00 

9.8000E-02 

».  rwr> ^  +  01 

4 .2.0  7  6  +  00 

6 . 4566+00 

9.9000E-02 

v  0  i 

4 .0676  +  00 

6.3316+00 

1 . OOOOF-O 1 

.  -i  •  6  f  k  c  f  0 1 

3.9806+00 

6.  194  6  +  00 
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Computer  Solution 

0~ a  60,000  psi 

STEP  =  5.02 

Steel 

TIME 

SUM 

V  ( 2 ) 

V  (  3  1 

0. 986 E- 06 

1 .0000E-03 

9.481E+0I 

1 .3  121.-02 

2 . 0000E-03 

9.470E+01 

0 . 23  3  E- 02 

9.06  If-- DO 

3.0000E-03 

9.4  69 1-  +  0 1 

1  .  1  7  2E-0  1 

4 .82 7E-04 

4 .  0000E-03 

9.462E+01 

2 . U70E-0  1 

1  .5  1  9 1:  -  0  3 

5. 0000E-03 

9.456E+01 

3.20  7 1 - 0  1 

3.6  /  1E-0  3 

6. 0000 E- 03 

9.449E+01 

4.0  7.3L-0 1 

7 .6081—0  3 

7. 0000 E- 03 

V  .  4  4  3  b  +  0 1 

6. I50E-0 1 

1  .  3991.-0  2 

6.0000E-03 

9.437E+0I 

7.924E-0  1 

2 .  .3681:— 02 

9 . 0000E-03 

9.4  30F  +  01 

9.B7  6  +  -0  1 

3.70  7E-02 

1 .0000E-02 

9.423E+01 

1 .  1 9  8  E +00 

0.66  71-02 

1  .  1  000E-02 

9 . 4  1  7  E  +  0 1 

1.42  3E+00 

M.  202E-02 

1.2000E-02 

9.4 10E+01 

1 . oOOE+OO 

1  .  14  7l:-0  1 

1 .3000E-02 

9. 404E+0I 

1 . 899F  +  00 

1 .0091  -0  1 

1  .4000E-02 

9.39/t+Ol 

2. 1 32E+00 

2.060E-0  1 

1 .5000E-02 

9. 390E+01 

2. 35SE+00 

2.6  7  7  E— 0  1 

1 .6000E-02 

9. 384E+01 

2.066E+00 

3.40  3E.-0  1 

1 .7000E-02 

i  9 . 3  7  7 1  +  0 1  i 

2 . 7601+00 

2. 9 021+ 00 

4 . 20  If  -0  1 

1 .8000E-02 

9.370E+01 

0 . 22  71—  0  1 

1 .  9000E-02 

9.363E+01 

3. 125E+00 

6.  3  36E-0  1 

2 . 0000E-O2 

V.356E+01 

3.286E+00 

7.080E-0  1 

2.  1  000E-02 

9.349E+01 

3.433E+00 

3.066E+00 

8.96  3 E-0  1 

2 . 2000E-02 

9 . 3  4  2  £  +  0 1 

1 .0491+00 

2.3000E-02 

9.3  308  +  01 

3. 686E+00 

1 .2  1  4  E  +00 

2 . 4000E-02 

9.328E+01 

3. 793F+00 

1 .  39  4  E  +  00 
1  ’.  086E+  00 

2.5000E-02 

9.32  1E  +  01 

3. 887E+00 

2. 6000E-02 

2 . 7000E-02 

9 . 3  1  4  tr  +  01 

3.9691+00 

1  '.  7921+00 

9. 30/E+01 

4.04 01  +  OO 

2  ’.009 1-  +  00 

2. 800 OE- 02 

9. 300E+01 

4 .  lOOfe  +  OO 

2  ’.  2  3  7  f  +  0  0 

2. 9000E-02 
3.0000E-02 

9.293E+01 

4.  I49E  +  00 

2 i  4  701  +00 

9.280E+01 

4  .  I9QE  +  00 

2;  72  31  +00 

3.  1  000E-02 

9.278E+01 

4 . 223E+00 

2  ’.9  781  +00 

3 . 2000E-02 

9.27  1E  +  01 

4 . 200E+GO 

3.24  1 1:  +  00 

3.3000E-02 

3. 40001-02 

9. 263 E +01 

4. 2711 +00 

3. 0081  +  00 

9.236F+01 

4.288E+00 

3'.  7801  +  0  0 

3. 5000E-02 

9.249E+01 

4 . 50  3E  +  00 

4-.O00E  +  00 

3. 6000E-02 

9.24  1E  +  01 

4.31  61: +00 

4  .  33 IE  +  DO 

3..  7000F-02 

9.234E+01 

4.329E+00 

4 . 608C+ 00 

3.8000E-02 

9.226E+01 

4 . 344  E  +  00 

4 . 88 21: +  00 

3 . 9000E-02 

9 . 2  1  9  h  +  0  1 

4.36 1E+00 

5. 10 4 E +00 

4.0000E-02 

9.2  1  1E  +  01 

4 . 383E+00 

0.4  2  )(-+  00 

4.  1  000F-02 

i  9.204E+01 

4.4 10E+00 

0.6  >131+00 

4. 2000E-02 

9.  196++01 

4 .4441+00 

5.7 3 /f .  + 00 

4,.  3000 E-  02 

9. 1 88E+01 

4 . 4 8SE+00 

6.  1841- +  00 

4..4000E-02 

9.  18  1E  +  01 

4 .03  0 1  +  0  0 

6 . 4  2  1  f-  »  00 

4. 5000E-02 

9. 1 73E+0I 

4. 090E+00 

6.648E+00 

4. 6000E-02 

9. 160E+01 

4. 660 E +00 

6.86  31  +00 

4. 7000 E- 02 

9. 157E+01 

4. 74 4 E +00 

7.06  7E  +00 

4. 8000 E- 02 

9. 149E+01 

4 .834F+ 00 

7. 2.0 Of  +  00 

<4 .  9000E-02 

9  .  1  4  2  fc  +  0 1 

4.934E+00 

7.  it  361  «  00 

5.0000E-02 

9.  1  34E  +  01 

5.046E+00 

7.6001  +00 
7  .  700 f  +00 

,5.  1 000E-02 

9.  1  26t  +  0 1 

0. I66E+O0 

5.2000E-02 

9.  1  18E  +  0I 

5.2968+00 

7.8867 +00 

5.3000E-02 

9.  1  1 OL  +01 

5. 4 34 E+ 00 

8. 00  71- +  00 

5. 4000F-02 

9. 102b +01 

0. OBOE +00 

Hi  i  ior +00 

5 . 5000F-02 

9. 094  E  +  0 1 

5.  7  3  3  C:  ♦  0  0 

H: 2091  +00 

5. 6000E-02 

9.085E+01 

5 . 89 1 E i 00 

8 . 2691-  +  OO 

■S.7000E-02 

9.077b  +01 

6.00  if i 00 

8.  30 61- +00 

5. 8000E-02 

9.069E+  01 

6.2 18E+00 

8.4  1  IF  +  00 

5.9000E-02 

9. 06  It +01 

6. 384 E  +00 

8 . 4  0  3E  +  00 

6.0000E-02 

9.053E+01 

6. 300E+00 

8.480F+00 

TIME 

6. 1 0006-02 

6.20006-02 
6. '>0001-02 

6. 40006- 02 

6. 50006- 02 

6.60006-02 

6.70006- 02 

6 . BOOOE-02 

6.90006- 02 
7.00006-02 
7. 1 0006-C2 
7 .200()E-02 

7.30006- 02 

7.40006- 02 

7.50006- 02 

7.60006-02 

7.70006- 02 

7.80006-02 

7.90006- 02 
6. 0000E-02 
8.  10006-02 

8.20006-02 

8 . 30006- 02 

8.40006- 02 
8.5000E-02 


SUM 

v  .  0 4 4  6  + 1)  1 
9.0  ■>(  l:  4-01 
>  .  0  2  7  >:  +  0 1 
9. 0  196  +  01 
9.0106+01 
9.0026+01 
8 .  9  ?  3  -  +  0 1 
6 .9846+01 
8.9766+ 01 
.  96  7l  +  01 
r,  .9586+01 
6. 9496+01 
8. 940++O1 
.  9  5  1  6  + 1)  1 
o. 9226+01 
8 .9  1  5  6  +  0 1 
3.9036+01 
•  .  8  94 1-  +  0 1 
8.385!  +01 
8.8756+01 
8 . 8666  +  0 1 
9.8566+01 
8. 8476+01 
8.8376  +  01 
8.82  76  +  01 


V(  2) 

6- 71 96+00 
6.8796+0" 
7  .  v  >  r’  >6+60 

7.  1756  +  00 
7.3146+00 
7.4436+00 
7. 9606+00 
7. 6636+ GO 
7. 7536+00 
7.8266+00 
7. 8836+60 
7.9236+00 
7.9446+00 
7.9476+00 
7.93 16+00 
7 . 8956+00 
7.84 16+00 
7.7696+C0 
7.6786+00 
7.8706+00 
7.44  56  +  00 
7.3056+00 
7.  130  6  +  00 
6.9826+00 
6. 8036+00 


8 . 6000E-02 

8.70006- 02 

8.80006-02 
8.9000F-02 
9.00006-02 
9.  10006-02 

9 . 20006- 02 

9. 30006-02 

9.40006-02 

9.50006-02 

9.60006-02 

9.70006- 02 

9. 80006-02 

9.90006-02 
1 .00006-0  1 


8.8175+01 
8.3086+0 1 
8.7986+01 
8.7886+0!  - 
8. 7786+0! 

<: .7686 +0! 

8. 7586+01 
3. 74  7+  +01 
3.7 576+01 
.72  76  +01 
.  7  1  76+01 
H.  706!  +01 
. 6966+01 
.6866 +01 
-.6  736  +  01 


6.6146+00 
6.4  16  6+00 
6.2126+00 
6.0 126+00 
5.8  1  86+00 
5.6526+0 0 
5.4536+00 
5.2826+00 
5.  1  19t  +  0() 
4.9656+00 
4.8206+00 
4.6836+00 
4.5546+00 
U. 4 34 6 +00 
4 . 3?  2 1:  +  00 


V(3) 

3.5056+00 
3.5166+00 
6.3  18  6  +  00 
8.5 1 26+00 
8.4^86+00 
8.4786+00 
3.4526+00 
8.4226+00 
58  76  +  00 
3.54 86 +00 
8.3076+00 
8 . 2646+00 
8.2  1  96  +  00 

8.  1  7  26  +  00 
8.  1246  +  00 
4.0746+00 
8. 0246+00 
7 .9756+00 
7.9?  1  6  +  00 
7.8676+00 
7.3126+00 
7 . 7556+00 
7.6956+00 
7.6326+00 
7.566F+00 


7.4956+00 
7.4206+00 
/.339E«-CC 
7.2326+00 
7.  1596  +  00 
7.0386+00 
6.9496+00 
6.8326+00 
6.7066+00 
6.5726+00 
6.4296+00 
6.2776+00 
6.  1  176  +  00 
5.9486+00 
5.7  705  +  00 
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APPENDIX  VI 


INSTRUMENTATION  SERIAL  AND  MODEL  NUMBERS 


Component 

Model 

Serial 

Velocity  meter  on  M. 

* 

L 

159 

Velocity  meter  on  M, 

* 

l 

222 

Velocity  meter  on  M, 

* 

3 

109 

Accelerometer  on  M2 

A5A-300-500 

4369 

Accelerometer  on 

A5A- 300-500 

4012 

Velocity  Meter  Control  122C 

Unit 

101 

Oscillograph 

5-119 

109025 

Amplifier  #1 

1-113B 

198BL10 

Amplifier  #2 

*  1-113B 

641BD09 

Power  Supply 

2-105 

245R4 

1000  cycle 

723-C 

2326 

*None  available 
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thesB24284 

Two-story  structure 


subjected  to 


002  01526  5 

KNOX  LIBRARY 


3  2768 
DUDLEY 


